arXiv:1508.04383vl [hep-ph] 18 Aug 2015 


Non-Simplified SUSY: r-Coannihilation at LHC and ILC 


M. Berggren 1 , A. Cakir 1,2 , D. Kriicker 1 , J. List 1 ,1.-A. Melzer-Pellmann 1 , B. Safarzadeh Samani 1 ’ 3 , C.Seitz 1 , 
S. Wayand 4 

X DESY, NotkestraBe 85, 22607 Hamburg, Germany 

Istanbul Technical University, Department of Physics Engineering, 34469, Istanbul, Turkey 
3 IPM, P.O. Box 19395-5531, Tehran, Iran 

4 KIT IEKP, CS Geb. 30.23 Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany 


DESY 15-145 


If new phenomena beyond the Standard Model will be discovered at the LHC, the properties of the 
new particles could be determined with data from the High-Luminosity LHC and from a future linear 
collider like the ILC. We discuss the possible interplay between measurements at the two accelerators in 
a concrete example, namely a full SUSY model which features a small ri-LSP mass difference. Various 
channels have been studied using the Snowmass 2013 combined LHC detector implementation in the 
Delphes simulation package, as well as simulations of the ILD detector concept from the Technical Design 
Report. We investigate both the LHC and ILC capabilities for discovery, separation and identification 
of various parts of the spectrum. While some parts would be discovered at the LHC, there is substantial 
room for further discoveries at the ILC. We finally highlight examples where the precise knowledge about 
the lower part of the mass spectrum which could be acquired at the ILC would enable a more in-depth 
analysis of the LHC data with respect to the heavier states. 


1 Introduction 


Supersymmetry (SUSY) |ljj6] is one of the best motivated theories beyond the Standard Model (SM). In 
the past, at LEP, Tevatron and the early LHC, results were mainly interpreted in the constrained minimal 
supersymmetric extension of the Standard Model (CMSSM) [7|[8] , while later LHC search results are mostly 
presented in a simplified model approach [9| |11| . Simplified models offer more flexibility for comparing to 
predictions of different theoretical models, though they might not describe reality well, as they usually 
contain only one decay with 100 % branching ratio. In contrast, a real SUSY signal might comprise a large 
spectrum of SUSY particles and the higher states of the spectrum may have many decay modes leading to 
potentially long decay chains. This means that the simplified approach in general does not apply beyond 
the direct production of the next-to-lightest SUSY particle (NLSP) and the interpretation of exclusion limits 
formulated in the simplified approach is non-trivial. Furthermore, also several production channels may be 
open, making SUSY the most serious background to itself. This becomes especially relevant for interpreting 
a future discovery of a non-SM signal. 


In this paper, we investigate how a well-motivated full-spectrum i?-parity 12 conserving SUSY model could 
lead to signals beyond the SM at the LHC and at a future linear collider like the ILC. A key feature of the 
investigated models is a r-NLSP with a small mass difference of about 10 GeV to the Xi , the lightest SUSY 
particle (LSP) of the model. Such scenarios were favoured by fits to all pre-LHC experimental data within 
the CMSSM [l3j, since the small mass difference allows to match the observed Dark Matter relic density 
via a sizeable r-coannihilation contribution. The T-coannihilation remains among the favoured regions of 
parameter space in the CMSSM and in the Non-Universal Higgs Mass Model (NUHM2) when including 
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LHC data 14 . Furthermore, light r’s will be an interesting candidate to enhance the rate for H —> 77 with 


respect to the SM expectation 15,16 


Within the context of constrained models, the masses of the fermion-partner particles (sfermions) in kine¬ 
matic reach of the ILC are excluded with high confidence. However, this exclusion is in most cases based 
on the strongly interacting sector, which in constrained models is coupled to the electroweak sector by 
GUT-scale mass unification. Without the restriction of mass unification, the part of the spectrum which 
is of interest to electroweak and flavor precision observables as well as dark matter, i.e. which is decisive 
for the fit outcome, is not at all in conflict with LHC results. This applies in particular to the t\ with a 
small mass difference to the LSP: Although first limits on direct r pair production from the LHC have been 


presented 17 , they rapidly loose sensitivity if the r is not degenerate with the e and and has a small mass 


difference to the LSP. In fact, the current limit on the r without other assumptions on the mass difference 
to the LSP than that it is larger than m T nor with any assumptions on the r mixing-angle comes from the 
DELPHI experiment at LEP and is My > 26.3 GeV 18 19 . Due the key feature of a small r-LSP mas 


difference and the resulting sizeable r-coannihilation contribution the series of CP-conserving model points 
considered here is called STC 


20 


Motivated by solving the naturalness problem within the general MSSM 21 , the SUSY partner particles of 


third-generation quarks in STC have been chosen to be lighter than those of the first- and second-generation 
squarks. When the first- and second-generation squarks and the gluino are rather heavy, > 2 TeV, the size 
of the total SUSY cross-section at the LHC strongly depends on the mass of the lightest top squark. We 
therefore consider in particular two model points, called STC 8 and STC10, whose physical spectra differ 
only by the mass parameter of the partners of the right-handed third-generation squarks at a scale of 1 TeV. 
The mass parameter of the right-handed third-generation squarks is set to 800 and 1000 GeV, resulting in 
physical masses of the top squark of mj i « 740 GeV and my ss 940 GeV, and of the bottom squark of 
mg « 800 GeV and mg ~ 1000 GeV, respectively. 

Figure [1J a) introduces the full mass spectrum of the benchmark scenario STC 8 , while Fig.Qb) zooms into 
the part of the spectrum accessible at the ILC with P cms = 500 GeV. The lightest higgs boson features 
SM-like branching ratios and has a mass in agreement with the LHC discovery within the typical theoretical 
uncertainty of ±3 GeV on MSSM higgs boson mass calculations. 



(a) (b) 


Figure 1: Full spectrum of STC 8 and decay modes with a branching fraction of at least 10% (a). The 
lower part of the spectrum of the STC scenarios, which features My x — M^o « 10 GeV (b). 


The dashed lines in Fig. [TJa) indicate those decay chains of the various sparticles which have branching 
fractions of at least 10%. The grey-scale of the lines indicates the size of the branching ratio. Only very 
few particles, namely the first and second generation squarks, the sneutrinos and the lighter set of charged 
sleptons have decay modes with 100 % branching ratio. 
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In particular the top and bottom squarks, but also the the superpartners of the uncoloured bosons, called 
electroweakinos in the following, have various decay modes, none of them with a branching ratio larger than 
50 %, but many with less than 10 %. This plethora of decay modes makes it challenging to separate the 
various production modes and identify each sparticle. 

The final goals of this study comprise the following questions for both LHC and the ILC in the example of 
the STC scenarios: 

• Which signature will lead to the first discovery of a discrepancy from the SM? How much integrated 
luminosity and operation time will be needed? 

• Which other signatures will be observable? 

• Which production modes of which sparticles contribute to this signal? 

• Can we tell how many sparticles are involved? 

• Which observables (masses, BRs, cross sections) can be measured and with which precision? 

• Can we show that it is SUSY? 

• Can the x® be identified as Dark Matter particle? 

In the next section, we will discuss the benchmark model and its phenomenology in more detail. In Section[3] 
we will review its possible discovery in different analyses at LHC and summarise the obtained simulation 
results, while we will do the same in Section [4] for the ILC case. In Section [5j we employ the combined 
results from LHC and ILC in order to investigate the questions raised above, before concluding in Section [6] 


2 Collider phenomenology of the STC scenarios 


In this section we present the parameters of the benchmark models in more detail, summarising the masses 
of the most important SUSY particles, their production cross sections and branching ratios. Based on this 
information, the phenomenology at the LHC and ILC will be discussed. 


A corner-stone of SUSY is that the couplings of the standard-model particles and their supersymmetric 
partners are the same, so that on tree-level the production cross sections only depend on the masses and 
mixing angles of the produced and exchanged sparticles. 


In the case of the LHC studies, we assumed E cins = 14TeV and an integrated luminosity up to 300 fb -1 . We 
also consider the high-luminosity upgrade of the LHC, the HL-LHC, assumed also to be running at U cms = 
14TeV, but delivering a total integrated luminosity of 3 ab -1 . 


For the ILC studies, the conditions presented in the Technical Design Report (TDR) [22] were used. This 
means that the bulk of the data would be recorded at E cms = 500 GeV, with an integrated luminosity of 
250 fb -1 per year. We extrapolate our results to the running scenarios recently published by the Joint 
Working Group on ILC Beam Parameters 23 . Since the ILC beam energy is tunable, we also consider the 
option of running at different, lower energies, and performing energy scans around thresholds. Very important 
are the opportunities offered by the electron and positron beam polarisation, for which the baseline design 
of the ILC foresees absolute values of 80% and 30%, respectively. As SUSY is a chiral theory, the possibility 
to have polarised initial conditions is a very powerful tool to disentangle different states, to enhance signal 
while reducing SM background, and to study helicity-dependent predictions of the theory. 


2.1 Mass spectrum and decay modes 

The masses of the most important sparticles in the benchmark points STC8 and STC10 are listed in Table[l] 
The two models differ mainly in the masses of the lighter bottom squark (bi) and the lighter top squark 
(ti). These sparticles are important for analyses at the LHC, where both would be accessible at the chosen 


3 



masses and will be searched for with dedicated analyses. A bottom squark with a mass as in the STC8 
model is on the edge of already being detected at the LHC with a centre-of-mass energy of 8 TeV and an 


integrated luminosity of 20 fb 


24 


25 


provided that the branching ratio of the direct decay b —► b%i is 


100 %. However, this is usually not the case for a full-spectrum model like the one we investigate here, unless 
the b is the NLSP. In particular, in our scenarios, the branching fraction of this decay is about 60 % in STC8 
and about 50% in STC10, as listed in Table [ 2 ] The analogous decay of the t has a branching fraction of 
only about 10 %. This will make it very hard to select a t sample in a specific decay mode as a prerequisite 
to identify a kinematic edge. 


Table 1: Sparticle masses for the models STC8 and STC10. The mass of the first two generation squarks (q) 
varies by a few GeV, their average mass is listed. We use £ when we refer to the first and second generation 
sleptons, while the third generation is listed separately. 


Sparticle 

Mass / GeV 

Sparticle 

Mass / GeV 


STC8 

STC10 


STC8 

STC10 

bi 

795 

1008 

n 

107 

107 

b 2 

1500 

1500 

r 2 

219 

219 

tl 

736 

944 

4 

196 

196 

t2 

1532 

1537 

X? 

96 

96 

g 

2042 

2042 

$ 

206 

206 

q 

2028 

2028 

$ 

410 

410 

4 

212 

212 

X ° 4 

425 

426 

4 

131 

131 

xt 

206 

206 

4 

198 

198 

xt 

426 

427 


Table 2: Branching ratios (BR) of the gluinos and the third generation squarks in the models STC8 and 
STC10. Only branching ratios above 1 % are shown. 


Decay 

BR 

STC8 

/% 

STC10 

Decay 

BR 

STC8 

/% 

STC10 

g -t tti 

37.8 

35.5 

g-t bbi 

39.1 

36.9 

g -t tt 2 

11.5 

13.8 

g -t bb 2 

11.5 

13.7 

4 -t tx? 

13.2 

9.9 

bi -> by? 

58.2 

51.4 

ti -t t*§ 

4.5 

4.3 

bi -> bx° 

3.1 

3.0 

4 -t 

22.4 

24.3 

bi -t byg 

10.7 

12.1 

4 ~11%4 

12.0 

15.1 

bi -> by° 

9.2 

10.5 

4 -> bxt 

10.8 

9.8 

bi -> txi 

5.1 

5.2 

4 -t b xt 

37.1 

36.5 

bi -> tx2 

13.7 

17.9 


The actual masses of the 4 and 4 in STCx have been excluded by the ATLAS experiment for the case of 
100 % branching ratio for the direct decay to the corresponding lepton and the LSP 26 . As can be seen in 
Table [s| the STCx branching ratios for these decays are quite close to 100%, thus this particular part of the 
spectrum is most likely excluded. It should be noted, however, that this large BR is a special case due to the 
small mass difference of only 6 GeV between the 4 and the xt / X 2 > which leads to a strong phase space 
suppression for the cascade decays 4 — > lx 2 and 4 —t viX 1 ■ With increasing mass difference, the branching 
ratio for the direct decay 4 —► lXi drops rapidly in favour of these cascade decays. Therefore, the results of 
this study remain highly relevant in the broader picture. Cases where the specific mass and branching ratio 
combination of the 4 plays a role will be pointed out. 


The STCx masses of the right-handed sleptons are by far not excluded, although their BR to lepton and 
LSP is 100%. This has to be attributed to their smaller mass difference to the LSP, which leads to softer 
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leptons and thus a significantly smaller acceptance. Among the other slepton decays listed in Table [3j it is 
interesting to note that the v T features 5 % of visible decays to TyW. 

Table 3: Branching ratios (BR) of the sleptons in the models STC8 and STC10. We use £ when we refer 
to the first and second generation sleptons, while the third generation is listed separately. 


Decay 

BR / % 

Decay 

BR / % 


STC8 

STC10 


STC8 

STC10 

£r —*■ lx 1 

100 

100 

4 -t rx°! 

100 

100 

4 -»• lx°i 

95.3 

95.4 

4 -t tx 1 

81.3 

81.5 

4 -»• lx°2 

1.7 

1.6 

4 TX° 

3.4 

3.4 

4 ->■ vixt 

3.0 

3.0 

4 -t v T xf 

6.2 

6.2 




r 2 T\Z 

9.1 

9.0 

4 -> ^Xi 

100 

100 

Vt -> v T x 1 

94.2 

94.6 




VT — > t{W 

5.8 

5.4 


Table 4 : Branching ratios (BR) of the charginos and neutralinos in the models STC8 and STC10. 


Decay 

BR / % 

Decay 

BR / % 


STC8 

STC10 


STC8 

STC10 

xt V- r*A 

0.2 

0.2 

X 2 -> efe* 

x ° ->• 

$ ->• 4^ 

2.1 

2.1 

Xi" -> w 

67.9 

67.5 

2.1 

2.2 

xt ->■ u e£ + 

6.6 

6.8 

73.2 

72.9 

xt ->■ 

6.6 

6.8 

X 2 -t v e v e 

5.7 

5.8 

xt -t WT+ 

11.3 

11.5 

X 2 -t Wn 

5.7 

5.8 

xf X?W+ 

7.2 

7.0 

X.2 -t V T V T 

9.5 

9.7 




X° -t x?z 

1.2 

1.1 

X 2 -+ e£i/ e 

4.6 

4.6 

->• x ? w * 

58.3 

58.4 

X 2 Ml 4* 

4.6 

4.6 

X° -t x?z 

10.3 

10.2 

X 2 T 2 ^r 

5.1 

5.1 

X° ->• X°Z 

23.2 

23.2 

X 2 -t *4e + 

1.7 

1.7 

x° x ? h ° 

2.2 

2.2 

X2 ->■ 

1.7 

1.7 

X°3 -t X° h ° 

1.2 

1.2 

X 2 -t VtT + 

2.5 

2.5 

Xa -t ' r 2 =r=F 

3.2 

3.2 

xt -> x ? w + 

7.8 

7.7 

X 4 *4*4 

4.3 

4.3 

X2 X°W+ 

28.3 

28.3 

X° 4*4 

4.3 

4.3 

X2 ->■ xfz 

25.0 

25.0 

X 4 -t v T v T 

4.3 

4.3 

X 2 -»• X?h° 

18.8 

18.8 

x° 4 -»• xfw ^ 

51.9 

52.0 




X '4 -> x?z 

2.3 

2.2 




X'4 ->• X°Z 

2.0 

2.0 




x°4 -> X ? h ° 

6.7 

6.7 




X'4 ->• X'2L 0 

15.8 

15.9 


Table [4] lists the decay modes of the electroweakinos. The largest branching ratios of about 70 % occur for 
t t\v t and xt 4 t > both always followed by t\ — ► t\ ?■ Since the T\ NLSP is only about 10 GeV more 

massive than the Xi LSP, the r leptons from the ri are typically very soft, making these decays challenging 
to detect. Although the mass splitting between x 2 (Xi ) and LSP is large enough to allow decays to on-shell 
Z (W*) bosons, these decay modes occur only at the level of a few percent due to the presence of the light 
sleptons. The heavier neutralinos Xi an d X4, however, feature sizeable branching ratios to xfW ± of about 
60% and 50%, respectively, while the x 2 decays most frequently to x^W*, xfZ and Xi h° with branching 
fractions of about 30%, 25% and 20%, respectively. Thus, although the t\ is the NLSP, in particular the 
heavier electroweakinos have sizeable branching fractions to other final states than the notoriously difficult 
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r-lepton. This also means that signatures with electrons or muons in the final state can originate either from 
slepton or electroweakino production. 


2.2 Production cross-sections at the LHC 


The cross sections for the main production processes of STC8 at the LHC are given in Table 0 tor E cms — 
14TeV, while Table [6] lists the cross sections which are considerably different in STC10. All cross sections 
have been calculated with Prospino2 27 28 . The particles produced with the largest production cross 
sections are the neutralinos and charginos. Nevertheless, their observation at the LHC is challenging due 
to their various decay chains via the ti, which lead to very soft r leptons due to the small mass difference 
between the f\ and the LSP. The cross sections for direct selectron and srnuon pair production is sizeable as 
well, such that these could be discovered in final states with two leptons and missing transverse energy. 


Table 5: Production cross sections for the benchmark model STC8 at LHC with E cms = 14TeV. The 
leading order (LO), next-to-leading order (NLO) and the Jv-factor between LO and NLO are shown. 


Process 

LO (fb) 

NLO (fb) 

-^NLO/LO 

PP -t g g 

0.18 

0.67 

3.58 

pp->qq 

5.1 

6 

1.16 

pp q g 

3 

5.4 

1.80 

pp->qq 

16.4 

25.4 

1.54 

pp bibi 

22.8 

38.3 

1.67 

pp -> b 2 b 2 

0.2 

0.37 

1.90 

pp -)• titi 

37.7 

62.7 

1.66 

pp — ► t 2 t 2 

0.16 

0.31 

1.94 

pp ->• 44 

15.5 

19.4 

1.25 

pp -)• 44 

32.4 

41.8 

1.29 

pp -)• 44 

19.4 

24.6 

1.26 

pp £4 

63.3 

79.6 

1.26 

pp -t 44 

59.1 

77.4 

1.30 

pp ->■ 44 

11.7 

14.6 

1.24 

pp 44 

34.7 

44.7 

1.28 

pp 44 

20.4 

25.9 

1.26 

pp —> TV t 

73.9 

93.6 

1.26 

PP -t XiXi 

0.34 

0.42 

1.25 

PP -t XiXj (except x°Xi) 

15.8 

19.6 

1.24 

pp -t xtxt 

585 

747 

1.28 

PP -t XkXm (except xtxf) 

17.3 

20.2 

1.17 

PP -t XiXt 

10.0 

12.9 

1.29 

PP -t X 1 X 2 

1.07 

1.36 

1.27 

PP -> X2X1 

1170 

1492 

1.28 

PP -t X2X2 

3.51 

4.33 

1.23 

PP -t X3X1 

6.10 

7.63 

1.25 

PP X3X2 

22.0 

27.0 

1.25 

PP -t XaXi 

3.45 

4.25 

1.25 

PP -t X4X2 

24.4 

29.8 

1.25 


In addition, in the direct decay bi —> by? leads to a characteristic edge in the contransverse mass distribu- 
, which depends on the masses of b and Xi- This edge is of high interest for determining the 


tion 


29 


30 


parameters of the model and could be observed if the direct decay mode can be selected with sufficiently high 
purity. In STC8, the cross section for bi pair production is sizeable, which gives hope to be able to observe 
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Table 6: Cross sections of the processes in the model STC10 at LHC with E cms = 14TeV. Only cross 
sections significantly different from those of STC8 given in Table [ 5 ] are listed. The leading order (LO), 
next-to-leading order (NLO) and the AT-factor between LO and NLO are shown. 


Process 

LO (fb) 

NLO (fb) 

A'nlo/lo 

pp ->• bibi 

4.4 

7.7 

1.74 

pp -S> titi 

7.08 

12.0 

1.72 


the edge. Due to the larger bi mass in STC10, the edge position is in a region with lower SM backgrounds, 
but on the other hand the cross section is lower. Thus, observing the edge will be challenging. Due to the 
rapidly decreasing cross section, we expect that bottom-squarks with masses well beyond 1 TeV cannot be 
seen at the LHC. 

The top squark with masses as chosen here will also be produced at the LHC, but it will be difficult to 
distinguish direct top-squark production from its production in the decay of the gluino (g). The masses of 
the heavier coloured sparticles, the gluino and the squarks of the first and second generation (q), are chosen 
such that they will also be produced at reasonably rates at the LHC, and dedicated analyses will be able to 
detect them with the full luminosity delivered by the LHC. 

The by far largest cross section in electroweakino production, above 1 pb, is obtained for X 2 X 1 production. 
This channel will therefore certainly be discovered at the LHC, by a multilepton search. The cross sections 
for other electroweak processes is lower, the largest cross section for neutralino-neutralino production appears 
for X 3 X 4 production with almost 14 fb. While the lighter electroweakinos could be well discovered at the 
ILC, the LHC searches would profit from the exact mass and cross section information of these in order to 
specifically search for the heavier electroweakinos that would not be accessible at the ILC at a centre-of-mass 
energy of 500 GeV. One example is the production of x 2 > which will be discussed later. 


2.3 Production cross-sections at the ILC 

The key feature of the STCx models for the ILC is the mass spectrum of the sleptons and the lighter 
electroweakinos, and thus at tree-level STC 8 and STC10 are identical from the ILC point-of-view. Figure [2] 
shows the polarised cross sectionf]for various STC processes in e + e“ collisions as a function of the centre-of- 
mass energy. In the part kinematically accessible at the ILC, they do not differ among the two models. While 
a few processes open up already below E cms = 250 GeV and thus would be accessible even when running near 
the higgs-strahlung threshold, a plethora of thresholds of slepton, sneutrino and electroweakino production 
appears between E cms = 250 GeV and 500 GeV. In most cases, these can be observed and even distinguished 
from each other in the clean ILC environment. The ability to operate at any desired centre-of-mass energy 
between 200 and 500 GeV (or even 1 TeV) and to switch the sign of the beam polarisations are unique tools 
to identify each of these processes. The low SM background levels allow in many cases a full and unique 
kinematic reconstruction of cascade decays. 

In particular at the ILC running at E cms = 500 GeV, all sleptons and the lighter set of electroweakinos 
of the STCx scenarios can be produced. x° and x° become accessible in associated production around 
-Ecms = 600 GeV. Pair production of x 2 appears at around E cms = 850 GeV. At this energy, also pair- 
production of X 3 and x° is possible; however because these two states are mainly higgsino, the rate is very 
low. At E cms = 500 GeV, the cross sections are sizeable, as can be seen in Table [7] Only one of the 
kinematically allowed processes, e + e _ —> t\T 2 , would have a production cross section below 10 fb for both 
beam polarisation configurations. The total SUSY cross section is well over 1 pb in both cases. 


1 Here, and in the following, 'P p + denotes the beam-polarisation configuration V{e . e~) = (p ,p + ), with p and p + 
given in percent (for brevity fully left (right) handed beams are denoted by L(R)). 
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E cms t GeV ] 



E cms [ QeV ] 


(a) (b) 

Figure 2: STC8 cross sections for sparticle production as a function of E cms at the ILC, separated in (a) 
sfermion production, and (b) electroweakino production. 


Table 7: Production cross sections for the benchmark model STC8 at the ILC, for different degrees of beam- 
polarisation. The ILC TDR beam-spectrum is used with a nominal centre-of-mass energy is 500 GeV. All 
channels accessible at this energy are shown. Channels with no detectable final states are marked with (*). 
In addition, the cross section for e + e~ -+ (e^e^) is 335.85 fb for Vr^r{Vl,l)\ for all other processes 

the cross-section vanishes for both Vl,l and Vr^r. 


Process 

Vr, l (fb) 

V l ,r (fb) 

E -(-80,-30 (fb) 

80,+30 (fb) 

e+e" -+ X?Xi(*) 

1203.57 

34.47 

705.29 

62.29 

e + e" -)• XiXi 

0.38 

259.37 

9.30 

151.74 

e+e" -> X 1 X 2 

11.97 

206.67 

14.24 

121.32 

e + e~ -+ X 2 X 2 

0.03 

115.71 

4.07 

67.69 

e + e~ -+ e L e L 

8.01 

84.07 

7.63 

49.46 

e + e“ -+ e R e R 

1313.73 

52.32 

770.36 

76.59 

e+e“ -+ /I L pt L 

8.05 

39.11 

6.08 

23.16 

e+e“ ->■ Jir/Ir 

222.47 

52.23 

131.97 

38.34 

e + e~ —> v T v T 

16.93 

22.02 

10.67 

13.47 

e + e~ v e v e {*) 

15.86 

973.97 

43.37 

570.33 

e + e~ -> 

15.93 

20.71 

10.04 

12.67 

e + e~ — > T\T\ 

244.33 

77.61 

145.65 

53.95 

e + e~ —> T 1 T 2 

5.46 

7.10 

3.44 

4.34 

e + e“ — ¥ t 2 t 2 

7.14 

26.10 

5.09 

15.52 

e+e“ -)• e L e R 

0.00 

0.00 

127.62 

127.62 


3 LHC projections 


The searches for new physics beyond the Standard Model at the LHC are either kept as inclusive as possible, 
or tailored to search for a specific scenario of new physics. We present here a representative selection of typical 


































studies, starting with search for and excess of large hadronic activity caused by the heavy new particles, in 
connection with large missing transverse energy due to the escaping lightest SUSY particles (LSPs). Such 
searches have been performed by both the ATLAS and CMS Collaborations based on data taken at 7 and 


8 TeV with and without b-tagging requirements 31 44 in the full-hadronic final state. Exclusive searches 


for top and bottom squarks rely on the use of special variables, as the production cross section for third- 
generation sparticles is about an order of magnitude lower than the one for first- or second-generation squarks 
of similar mass and the signal is hidden below a large background by SM top-quark production. We discuss 
here a full-hadronic search for direct bottom squark production, similar to searches performed by both the 
ATLAS and CMS Collaborations based on data taken at 7 and 8 TeV 24 25 45 , where it is assumed that 


the pair-produced bottom squark decays directly to a bottom quark and the lightest neutralino which is 
the LSP. Furthermore, we present the results of search for direct top squark production in the single-lepton 
channel, also similar to previously published analyses by CMS and ATLAS 46-48 . Multilcpton searches 


are sensitive to decays of the electroweak-produced sparticles, e.g. in case of neutralino-chargino production, 
as discussed for 7 and 8 TeV by the ATLAS and CMS collaborations. The recent results are interpreted in 
the SMS approach which reflects a best case scenario 49-51 , assuming a branching ratio of 100% for one 
specific decay. In nature the 100 % are not realised, as also shown for our benchmark models in Table [2] 
For such cases the bounds for the chargino and neutralinos are weaker. We also present the results of a 
multilepton search and the possible interpretation in case of a signal by the discussed model. 


The detector response is simulated using the Delphes 3.0.9 fast simulation program 52 , both for signal and 


background events. This Delphes version has been used for the Snowmass studies and is able to include 
pileup from minimum bias collisions that are randomly selected from a file containing inelastic proton-proton 
interactions produced with Pythia6 53 . These events are randomly distributed along the beam axis (also 


called 2-axis) according to a Gaussian distribution with a width of 5 cm. If the 2-position of a pileup vertex 
is less than the 1 mm from the primary vertex (corresponding to the resolution), the pileup interaction is 
not separated from the primary vertex, and all particles from both the pileup and primary interactions are 
included in the object reconstruction. For pileup interactions with a larger 2-vertex difference to the primary 
vertex, the subtraction of charged pileup particles within the tracker volume is applied with an efficiency of 


unity. The Fast Jet area method 54 is applied to correct measurements of jets and energy in the calorimeters 


for the contribution from neutral pileup particles and charged pileup particles outside the tracker acceptance. 


About 10 to 100 million events per background process that were produced for the Snowmass effort 
with Madgraph5 


55 


56 , including up to four extra partons from initial and final state radiation, matched 


to Pythia6 for fragmentation and hadronisation, are used in this paper. The background cross section is 
normalised to next-to-leading order (NLO) in the background production process, which is based on the work 
in preparation for the Snowmass summer study 2013 and discussed in more detail in Ref. 57 59 . While 
we studied all the major sources of background events, background processes with low cross sections that 


might become relevant at 3000 fb -1 are not included. The signal samples are generated with Pythia6 and 
passed through the Delphes simulation. For Pythia6 the tune Z2* 
are calculated at NLO with Prospino2 


27 28 


60 is used. The signal cross sections 


Assuming systematic uncertainties of the same order as in the existing 8 TeV analyses, we determine for 
each search the discovery sensitivity, using the Binomial significance Zb i 61 63 in Roostats 64 . Here, the 


sensitivity is calculated in a frequentist way in one-sided Gaussian standard deviations, performing a hy¬ 
pothesis test between background-only and signal-plus-background, where the uncertainty on the background 
estimate is taken as Poisson distributed. 


3.1 Full-hadronic search 

Heavy squark and gluino production in R-parity conserving SUSY scenarios can lead to long decay chains 
with multiple jets and therefore a large amount of hadronic energy, and large missing transverse momentum. 
A typical search for such a scenario is based on the variable H t, the scalar sum of the momenta of all 
jets with pt > 50GeV and |7y| < 2.5, and missing hadronic transverse energy (MHT), which is defined as 
absolute value of the negative vectorial sum of all jets with px > 30 GeV and |7y| < 5. The SM background to 
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this SUSY search arises mainly from the following processes: Z(z^) + jets events, and W (lv) + jets events 
from W, or tt + jets, where at least one W boson decays leptonically. The W(lu) + jets events pass the 
search selection when the e/p escapes detection or when a top decays hadronically. QCD multijet events also 
contribute to the background when jet-energy mismeasurements or leptonic decays of heavy-flavour hadrons 
inside jets produce large MHT. However, the QCD background generally becomes negligible at very high 
MHT as required here. 


The analysis follows the baseline selections motivated by the 8TeV analysis 41 , requiring at least three 


jets with px > 50GeV and \rj\ < 2.5, H x > 1000 GeV and MHT > 500 GeV. In order to remove the QCD 
background, events are required to satisfy the following cuts for the azimuthal angle difference between the 
leading jets and MHT direction: |A< p(j n , MHT)| > 0.5 for n = 1,2 and |A</(j 3 , MHT)| > 0.3. Events are 
vetoed when they contain isolated muons satisfying px > 10 GeV and |?y| < 2.4 or electrons with px > 10 GeV 
and |?y| < 2.5, which suppresses mainly the tt and W(li/) + jets background. The baseline cut flow is given 
in the Appendix, in Table [9j 


Table 8: LHC all-hadronic inclusive search: Background and signal event yields corresponding to 300 fb 1 
The notation ”V” refers to W, Z and 7. Four signal regions (SR) are shown here. 


Signal regions 

tt +jets 

V+jets 

VV+jets 

Top+jets 

Total SM 

STC8 

STC10 

SR A: 

H t > 2000 GeV; 
MHT > 1500 GeV 

0.29 

43 

2.1 

0.01 

45 

45 

35 

SR B: 

H t > 2000 GeV; 
MHT > 1000 GeV; 
n(Bjets) > 2 

4.4 

14 

0.92 

0.15 

20 

51 

35 

SR C: 

H t > 3500 GeV; 
MHT > 1000 GeV 

0.49 

9.5 

0.54 

0.01 

10 

8.4 

7.6 

SR D: 

H t > 3000 GeV; 
MHT > 1500 GeV 

0.1 

6.3 

0.41 

0.002 

6.8 

8.2 

7.2 


The recent searches by the CMS Collaboration split the events from the baseline selection into several 
exclusive search regions according to their Ht, MHT and b-tag multiplicity, while for this study we keep 
only four well-motivated signal regions listed in Table [8] The most promising signal region targeting the 
inclusive production of heavy gluinos and squarks of all generations is SR A, where we require H x > 2000 GeV 
and MHT > 1500 GeV. Figure [3j a) and (b) shows the Ht and MHT distribution for this selection. Two 
other signal regions, SR C and SR D, are characterised by higher ifx requirements and reject too much 
signal, which could be taken as a hint on the (not too high) squark or gluino mass. In SR A we find 45 (35) 
signal events for STC8 (STC10) over 45 background events for an integrated luminosity of 300 fb -1 . About 
80 % of these signal events are first- and second-generation squarks. Assuming a systematic uncertainty of 
20 %, we determine a discovery sensitivity of about 3 cr. 

Another signal region, SR B, defined by Ht > 2000 GeV, MHT > 1000 GeV and TVb-tags > 2, is tailored to 
SUSY signals with a light third generation, being sensitive to either gluinos decaying through (virtual) top 
or bottom squarks or to directly produced third-generation squarks. The results for this signal region are 
shown in Fig. [3](c) and (d). We find 51 (35) signal events in STC8 (STC10), expecting 20 background events 
for 300 fb -1 . A discovery sensitivity of 5 a will be reached with 200 fb -1 , assuming a systematic uncertainty 
of 20%. The signal consists mainly of gluino-associated (g g and g q) and direct heavy-squark production 
(about 65% of the final events). 

Further studies of the kinematic variables may shed more light on the nature of the new physics seen in this 
scenario and are discussed in the following sections. 
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Figure 3: LHC all-hadronic inclusive search: Ht and MHT for signals and SM backgrounds after all 
selection requirements are applied. (a),(b): inclusive search region SR A; (c), (d) signal region SR B 
including the requirement of at least two b-tagged jets. 


3.2 Search for direct bottom squark production in the final states with two 
b-quark and missing energy 


In this section, we investigate the discovery potential for third-generation squarks in the final jstate with 
two b-jets and missing energy in the LHC. Of particular interest for this search is the decay bi —» by^, 
which is the dominant decay mode of bottom squark in the investigated models, with branching fractions of 
more than 50 %. Assuming that the bottom squarks are pair-produced, a final state containing exactly two 
b-quarks and two neutralinos is expected for about 25 % of the signal events. As the masses of the xt and 
the Xi are not degenerate, the contribution of the top squark to this final state is small the STCx models. 

Two jets which originate directly from bottom squark production should exhibit sizeable transverse momen¬ 
tum. Figure |4(a)| and (b) display the px of the two leading jets for all the standard model processes and the 
signal. When the bottom squark decays directly to a b quark and a neutralino, the resulting b-jet px is on 
average harder than for SM and other SUSY processes. Therefore, we require the events to contain exactly 
two central jets with \g\ < 2.4 and px > 300 and 200 GeV, respectively. Events including any additional jet 
with px > 70 GeV are rejected. In addition, both leading jets are required to be identified as originating 
from a bottom quark. 
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Figure 4 : LHC bottom squark search: (a) pt of the leading jet for events with exactly two central jets; 
(b) pt of the next-to-leading jet after applying the leading jet pt requirement; (c) E™ lss after requiring for 
exactly two b-tagged jets; (d) m t after requiring for exactly two b-tagged jets and E™ lss > 450 GeV. All 
signal events where at least one bottom squark do not decay to bi -4 by) 1 , plus all other SUSY processes 
are classified as “STC other”. 


By vetoing events with at least one lepton with pt > 10 GeV, we suppress the main standard model back¬ 
ground processes such as tt and W + jet s. The distribution of the missing energy for the events passing the 


above requirement is shown in Fig. 4(c) We require missing transverse energy to be E!f lss > 450 GeV. 


The tt and W + jets backgrounds are suppress further by requiring the minimum invariant transverse mass 
tot of one of the two leading jets (jeti) and the missing transverse momentum to be greater than 500 GeV. 
The invariant transverse mass is a quantity that is often used in LHC analyses, and it is defined as function 
of the transverse momenta of two objects a and b as follows: 


m T (pT,PT) = ^PtPt (! ^ cos(A</>(p!,f4))) . 


The distribution of min 


™ zumi 


s\ i tp : 

) > TTIt (Pt 5 -Ca 


'miss 

T 


kinematic edge at the mass of the top quark when the jet and E™ 1 
top quark. A outflow is shown in the Appendix, in Table |10| 


(i) 


is expected to have a 


originate from semileptonic decay of a 


as shown in Fig. 


4(d) 
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With an integrated luminosity of 300 fb” 1 , we find 74 (44) signal events in STC8 (STC10) for a background 
yield of 28 events. 73 % of the selected signal events in the STC8 scenario are from direct bottom squark 
production with the desired decay, bi —» by?. In STC10, this number reduces to 58% due to higher 
contamination from decays to b quarks and yij, yg, or Xi- In both scenarios, we observe about 5% SUSY 
background due to ti decays. The other bottom squark decay modes and all events originating from top 
squark decays are treated as SUSY background in this analysis. Assuming a systematic uncertainty of 15 %, 
this channel reaches a discovery sensitivity of 5er with about 90 and 300 fb -1 of data for STC8 and STC10, 
respectively. 


3.3 Search for direct top squark production in the single-lepton channel 


In this section, we discuss the search for direct top squark pair production. Previously conducted searches by 
the CMS collaboration 48 focus on simplified SUSY models, where only the process of interest is considered 
while all other sparticles masses are assumed to be out of reach. In these simplified models two top squarks 
are produced, which decay into either tx? or byy, with varying branching ratios. The two STC models 
considered here, however, only a very small fraction of events where two top squarks are produced exhibit 
the sought after decay process. As an example, in the STC8 model only 1.6% of tit* events are expected to 
have a decay mode tit J —► (t,Xi)(t,Xi)- 


Additionally, bottom squark pair production enters as a sizeable intrinsic SUSY background when one of the 
bottom squarks decays to by 1 / and the other one decays to ty x . Gluino and squark production also enter 
as an intrinsic background with top quarks in the decay chain resulting in similar signatures as expected for 
process of interest. Thus, we face two major challenges in this analysis, one being large SM background and 
the other one being background from other SUSY processes. 


The analysis method follows the aforementioned search performed by the CMS collaboration at 8 TeV 48 


but with tighter selection requirements. We require a single isolated electron or muon with px > 30 GeV 
and |?7| < 2.4. Events are vetoed if there are additional isolated leptons with px > 20 GeV. In addition, 
we require at least five jets with px > 40 GeV and \r]\ < 2.4, which enhances the fraction of titJ events 
with respect to bibj events. One or two of these jets must satisfy at least a medium b-tag requirement. To 
further reduce the SM background, we require E™ lss > 400 GeV. 


Additionally, we introduce an angular variable min A<p, the minimum azimuthal angle between the leading 
or sub-leading jet and the E™' ss . For this variable we require events to have a value greater than 0.8 in order 
to reduce backgrounds from SM processes. Another variable that aids in reducing backgrounds is centrality, 
defined as the sum of the px of the lepton and jets divided by their total momentum ^ jC je^ P (p)+Iepton(p^ ■ 
For SUSY events we expect this variable to be shifted towards higher values, while SM backgrounds are less 
central. Events are selected that satisfy centrality > 0.6. 


After requiring the transverse mass, mx, calculated with Equation[T]for the system consisting of the lepton px 
and the missing transverse momentum vector, to satisfy mx > 260 GeV, the background arises predominantly 
from two sources: tt events in which both W bosons decay leptonically but one lepton is lost, and diboson 
events. In order to suppress the tt background, we require m^X,, defined as the minimum “mother” particle 
mass compatible with all the transverse momenta and mass-shell constraints 65 , to be above 260 GeV. By 


construction, for the dilepton tt background without mismeasurement effects, has an endpoint at the 
top quark mass, while for semi-leptonic tt events and signal it has a large tail. Figure [s] shows the A</>, 
centrality, mx, and the rn^ 2 distributions after all previously mentioned selection criteria are applied, except 
on the variables themselves. A outflow can be found in the Appendix, in Table El 


After all selection requirements, we find for 300 fb” 1 155 (STC8) and 113 (STC10) signal events with a 
very low background of 28 events. Assuming a systematic uncertainty of 15%, the STC8 model could be 
discovered with 5cr with an integrated luminosity of 25 fb” 1 , and STC10 with 40 fb” 1 . While the search 
targets direct top squark pair production, the final composition of the selected events contains only about 42% 
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Figure 5: LHC top squark search: Comparison of (a) A <f>, (b) centrality, (c) m t, and (d) rn^ 2 for signals 
and SM backgrounds after all selection requirements are applied except on the variables themselves. The last 
bin in in the variable m^ 2 holds all events where the minimisation, which is performed during the calculation 
of the variable, does not converge within a reasonable amount of steps. 


of this type for STC8 and 23% for STC10. The additional events originate from squark-gluino production 
and bottom-squark pair production. 


3.4 Search in the multilepton channel 

In the STCx models, the cross section for direct %2 Xi production in proton-proton collisions amounts to 
more than 1 pb. The golden channel for the discovery of this process is the multilepton channel, with final 
states containing three or more prompt leptons. The leptons for the signal are produced by slepton mediated 
decays or by leptonic decays of W or Z bosons which are produced in the decay chain. However, in the STC8 
model the t) is the NLSP and almost mass degenerated with Xi> leading to soft leptons in the dominant 
decay modes of \2 an d xt decays via ti (c.f. Table [4j. Most of these leptons will not pass the selection 
criteria, and thus we expect only a small sensitivity to these decay modes. 

In the following we define lepton to be either an isolated muon or an isolated electron, which includes the 
leptonic decay modes of the r leptons. As hadronic decays of r leptons are not optimally modelled in 
the used Delphes version, we do not consider these here. All leptons must have a transverse momentum 
Pt > 10 GeV. In order to comply with the trigger, the leading (sub-leading) lepton must satisfy pt > 25 GeV 
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(15GeV). We present here the analysis of the three-lepton final state requiring exactly three leptons. The 
analysis of events with four or more leptons does not increase the sensitivity and is not discussed further. 

The WZ and ZZ production are summarised as VV background. Non-prompt backgrounds cover all events 
which include leptons from misidentified objects (also known as ’fake’ sources), or leptons which are not 
produced by the hard scattering. For example, dileptonic decays of tt, where one of the b-jets produces an 
isolated lepton, leads to final states with three leptons. All other processes which contribute to the three- 
or four-lepton final states, e.g. SM higgs-boson production or triple-boson production, are summarised as 
rare backgrounds. The signal is subdivided into four production mechanisms, the direct chargino-neutralino 
production of the second-lowest mass electroweakinos and the production of at least one higher-mass 

chargino, xf > summarised as X2 Xm- Other direct production modes of charginos and/or neutralinos are 
comprised in ’other EWK’, while we label all events from other sources than electroweakino production, such 
as slepton-pair production or production of coloured SUSY particles, as ’noEWK’. 

The main background in the three-lepton final state originates from WZ production, where both bosons 
decay into leptons. The transverse mass, calculated with Equation [l] for the missing transverse energy and 
the pt of the lepton from the W boson, is used to suppress this kind of background. The second most 
important background are three-lepton events from non-prompt sources. Most of those events stem either 
from tt or Drell-Yan events. The latter are suppressed by the requirement of large E7p lss , whereas the tt 
events are reduced by a b-jet veto. In order to maximise the sensitivity we use a three-dimensional binning 
in the variables m g+g~ (invariant mass closest to mz), U™ lss and tot- In Fig. |6j the expected distributions 
for SM and STC8 are shown as well as definition of our search regions. The events are further separated by 
the value of the invariant mass of the opposite sign same flavour lepton pair: m g+g- < 75 GeV, 75 GeV < 
m f.+£- < 105 GeV and m g+p- > 105 GeV. 

Overall we define 45 independent search regions, which are listed in the Appendix in Table [12] and which 
all contribute to the final sensitivity. The applied uncertainties include lepton uncertainties (3% per lep¬ 
ton), uncertainties on the E™ lss shape (15%-25%) and uncertainties of the MC statistics (0%-30%). When 
combining all search regions, we can discover STC8 with less than 200 fb. The most sensitive individual 
regions are at medium wit (200-400 GeV) and medium E™ lss (200-400 GeV), either requiring mg+g- to be 
around the Z boson mass, which is mainly driven by direct X 2 X 1 production followed by X2 — t ^Xi, or 
for m £4- g- higher than the Z boson mass. In the latter signal region long-chain decays of the higher-mass 
neutralinos and charginos give important contributions, whereas the role of decays of the X2 via smuon or 
selectron is reduced due to their kinematic edge at m g+g- < 110 GeV. 


4 ILC projections 


Due to the clean conditions at an e + e collider, the searches for SUSY at the ILC are exclusive, with separate 
analyses adapted to each specific channel searched for. 


Out of the plethora of exclusive processes, in particular pair production of sleptons has been analysed in 
context of this paper, assuming an integrated luminosity of 1000 fb -1 at E cms = 500 GeV, which corresponds 
to four years of data taking at the design luminosity of the ILC. The data is assumed to be collected in 
two equal-size samples, with beam polarisations of V+so ,-30 or "P-so.+sOj respectively. The beam-conditions 
described in the ILC TDR 22 are applied for the simulation. We use the same SM background samples as 
for the detector benchmarking, presented in the ILC TDR 66 . These samples are produced with v.1.95 of 
the Whizard 67 event-generator. For one background category, 77 processes, the statistics in the TDR 


sample was insufficient, and a dedicated simulation, using the same version of Whizard is done to extend 
also this sample to correspond to the above-mentioned integrated luminosity of 500fb _1 per polarisation. 
Under these conditions the total expected SM background is 1270 million events. The SUSY signal is also 
generated with Whizard 1.95, using the same settings as for the SM samples, with the exception that 
the study of the threshold scans of and eR, where the signal is generated using Pythia6.422 53 . The 
detector simulation used is the fast simulation program SGV 68 , adapted to the ILD detector concept at 
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Figure 6: LHC Multilepton search: E™ lss (a), mx (b) and m f+f- (c) distribution for three lepton events. 
Figure (d) illustrates the numbering scheme of the search regions. Detailed results for these are given in the 
Appendix, in Table 
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ILC 66 . In this program, the detector geometry is described in a complete but simplified model, and is 


used to calculate the helix parameters and covariances for all charged particles individually. The calculation 
includes description of the point-resolutions at each layer, the magnetic field, and the effect of multiple 
scattering. In addition, brems-strahlung and 7 conversions in the tracker material is simulated at each layer. 
The measurements by the calorimeter-system are in a first pass simulated parametrically for each particle 
reaching it. In a second pass, the individual showers are split or merged by a procedure adjusted to as well 
as possible emulate the performance of the particle-flow algorithm as implemented in the the full simulation 
and reconstruction of the ILD detector. 


Beyond slepton pair production, we summarise the results from previous studies performed on similar bench¬ 
mark models and discuss them in the context of STCx. Most of these studies have been performed in full 
GEANT4-based simulation of the ILD detector concept, using samples simulated for the ILD Lol 69 . 


4.1 Analysis of direct electroweakino production 

At the ILC, direct electroweakino production occurs at rates in the order of 100 pb, c.f. Fig. [2] At a centre-of- 
mass energy of 500 GeV, XiXi> X1X2, X2X2 and Xi Xi are kinematically accessible in the STCx scenarios. By 
measuring their polarised production cross sections and the masses of the involved sparticles, the parameters 
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of the electroweakino sector (Mi, M 2 , g, tan/3) can be determined 
experimental observations. 


70 , given sufficient precision of the 


As a particular challenge of the STCx models, the y? an d Xi decay dominantly via the ri NLSP. Nevertheless, 
almost all electroweakinos have branching fractions at the few percent level to other final states than the 
notoriously difficult r lepton. The use of beam polarisation and tunable E cms will further enhance the power 
of the observations and allow to disentangle production modes with very similar final states. 


4.1.1 X 1 X 2 production 


The dominating decay X 2 ~► DT, with a branching ratio of about 70 %, leads to the same final state content 
(t + t _ x?x?) as for 7) pair production. However, they can be disentangled by their different beam polarisation 


dependency and kinematic properties. The di-r invariant mass can be employed to measure the mass of the 
X 2 - This channel has been studied based on four-vector smearing, indicating that precisions of 1-2 GeV could 
be achievable 


71 


In our scenario, the di-r final state also receives background from X 1 X 1 ► t\v t t\v t , 
which features a similar cross section x BR at E cms = 500 GeV. The two processes can in principle be 
disentangled by the different kinematic features, but the achievable resolutions would need to be studied. 
However, thanks to the tunable centre-of-mass energy of the ILC, this problem can be avoided altogether by 
collecting data below the threshold for chargino pair production, in our case e.g. between E cms = 350 and 
400 GeV, or by scanning the X 1 X 2 production threshold. 


More recently, it has been shown in full GEANT4-based simulation of the ILD detector that the contribution 
from X 1 X 2 MMX 1X1 j which is about a factor of 10 smaller, leads to competitive mass and cross section 
measurements. Figure [7] shows for instance the invariant mass spectrum of the two muons before ((a), signal 
is scaled by factor 100) and after event selection (b). From this channel alone, the mass of the y 2 can be 
determined to a precision of about 1 GeV for an integrated luminosity of 500 ftV 1 collected with Vso,+eo, 
depending on the assumed precision for the mass of /Ir and y? 72 . 
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Figure 7: Determination of the y? ma ss from the di-muon invariant mass spectrum, with the full spectrum 
for inclusive di-muon sample (a) and a zoom into signal region after dedicated selection (b). From 


72 


Within the same study, the corresponding uncertainty on er(e + e _ —>• Xi,\ 2 ) x BR(x 2 * MMX?) bas been 

determined to about 20%, while the precision on a(e + e~ —>• X1X2) x BR(x 2 TT Xi) bas been estimated 
to 2 % 71 . We thus conclude that the corresponding studies of the di-r channel should be repeated with 


up-to-date simulation of the expected detector and accelerator performance. 
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4.1.2 XiXi production 


As discussed above, the final states of xt pair production will be dominated by di-r plus missing four- 
momentum final states from xt —> T\V r or Xi —t v t t followed by an invisible decay of v T: thus be similar to 
those from XiXti- While at E cms = 500 GeV the cross section x BR should be clearly measurable above the 


X 1 X 2 background, especially once the latter is known from running below the Xi Xi threshold, kinematic 
mass reconstruction in this channel needs further study, like in the case of XiX 2 - However, there are several 
alternatives to measure the mass: 

On one hand, the xt features a branching ratio of about 7% to x?Events where this decay, followed 


by W 


qq', 


occurs for one of the Xi ■, while the other could e.g. decay to t\v t , give a very unique signature. 


The edges in the energy spect rum of the W bosons can then be used to determine the mass of the xt if the 
Xi mass is known (see Section 4.2). This reconstruction method has been studied in full detector simulation 
by ILD 69 73 and SiD 74 for a SUSY scenario where the sleptons are heavier than the electroweakinos 
XiW and x .2 XiZ have branching ratios close to unity. These studies achieved 


10 3 Xi —•*■ XiW decays, assuming the Xi mass 


and thus both Xj 

mass resolutions of about 1.5% for the Xi based on 125 
is known from another source. In STCx, with an integrated luminosity of 500 fb -1 , only 11 • 10 3 decays 
would be available and the Xi mass would be known to permille level precision from slepton pair production, 
However, backgrounds would also be lower, since a) the branching ratio for xti 


c.f. Section 


4.2.1 


_ X°iZ 

is 100 times lower than in the original scenario of these studies and b) the dominating SM background 
does not originate from fully hadronic W pairs anymore, but from semi-leptonic ones, where the charge 
of the lepton in conjunction with the forward-backward asymmetry provides an additional, very effective 


suppression mechanism 75 . But even neglecting these expected benefits, the pure scaling according to the 
number of decays yields a projected uncertainty of 5% on the xt mass. With the full running program of 
the ILC 


23 , this would shrink to about 2.5%. 


In addition, the decays Xi _■ y GT- and Xi ~t give important information not only on the Xi mass, 
but also on the mass of the mostly invisible v. Since the mass difference between xt an d ^ is only 10 GeV, 
lower as well as the upper edge of the lepton energy spectrum are significantly below the lower edge from 
pair production of the left-handed sleptons (c.f. 14.2.11), at 8 and 16 GeV, respectively. Pair production of the 


right-handed selectrons is heavily suppressed by the appropriate choice of beam polarisation. Solectron and 
smuon backgrounds can be reduced further by selecting chargino decays to different lepton flavours in the 
two decay chains of an event. While t\ pair production leads to different flavour leptons in the final state, 
the decay leptons of the r decays will be distributed over a wide range of energies (c.f. Fig[l2j) , whereas the 
signal leptons lead to a sharp-edged, narrow box with the above edge positions. With the branching ratios 
of STCx, about 20-10 3 decays to Vfi will be available from 500 fhU 1 of data. Since this is of similar size than 


the available statistics for e.g. the hl (c.f. sec 4.2.1), with steeper edges due to the smaller range of lepton 
energies, a precision of 1% or better should be achievable for both the xt and the v mass. 


4.1.3 Radiative X 1 X 1 production 


Additional information can be gained from the analysis of radiative X 1 X 1 production, employing a similar 
technique as in WIMP dark matter searches in mono-phot on signatures 
GEANT4-based simulation of the ILD detector concept in 


77 


76 . This has been studied in full 


With an integrated luminosity of 500 fb 


equally split between all four beam helicity combinations for |P(e“,e + )| = (80%,30%), the unpolarised 
cross section as well as the cross section for purely right-handed electrons and left-handed positrons could 
be determined with a precision of about 5 % each. The mass of the x? could be determined to about 2 GeV, 
completely dominated by a rather conservative estimate of the systematic uncertainty due to the limited 
knowledge of the shape of the beam energy spectrum. 


As can be seen in Fig.[8j the statistical precision, here from a template fit to the energy spectrum of the ISR 
photons, is much lower already with this rather modest luminosity assumed in the study. A more detailed 
analysis of the impact of the beam energy spectrum should be performed in the future. 
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Figure 8: Determination of mjo from a template fit to the photon energy spectrum in e + e 
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4.2 Analysis of direct slepton production 


By studying slepton-pair production in the continuum (i.e. well above the threshold), one takes advantage of 
the known initial state at the ILC: In general, the endpoints of the energies of the standard model particles 
X and Y (both assumed to have negligible mass) in the proces^] e + e _ —► XY —> XYxiXi can be found to 
be 



Mr f (MfO \ 2 \(& ab + 1^1 \ 
2 [ v Mr) ) { M~ (-) M~ ) 


( 2 ) 


where i is either X or Y. Note that due to momentum conservation, p l ~ b = —p l ~ b . If the two sparticles 
have the same mass, so that M^ = My = M~, also E l ~ b = E l ~ b = E cms / 2, and \pj ab \ = J{E C ms/2) 2 — M~. 
Hence, in the case of slepton pair production, e + e“ 1X1, one finds that 



i.e. by determining these two endpoints, and using the knowledge of E cms , both M^o and Mj can be deter¬ 
mined. Even in the case where the slepton is a r, the endpoint of the spectrum of the r decay products can 
be used to determine My, provided M^o is known from e.g. the or and spectra, measured from the same 
data. 


The key characteristics of slepton production and decay, which single them out from the background, are: 

• only two leptons in the final state 

• large missing energy and momentum 

• large acollinearity, with little correlation to the energy of the lepton 

• central production 

• no forward-backward asymmetry 

Different backgrounds dominate for the lighter and the heavier sleptons: for the ri, the background from 
photon-photon processes is important, while WW —> Ivlv is less so; the opposite is true for the T2 , ?l and /7 l- 

2 Usually, X and Y are particle and anti-particle. 
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For e and ju, the lower endpoints of the spectra are high enough - both for the right-handed and left-handed 
states - that the photon-photon background does not pose a problem. The mass measurement relies on the 
detection of edges or endpoint in the spectrum of the SM decay products. Therefore, the SUSY background 
is less important, since it is only important that it is small and/or flat in a narrow region around the edges. 
In case for e and /t the SUSY background is indeed found to be both small and flat in these regions. In the 
case of the r, the SUSY background is dominated by X 1 X 1 and X 2 X 2 production with cascade decays over r 
sleptons, which all yield an upper kinematic limit of the produced r leptons well below those of both 7/ and 
r 2 pair production, so they have little influence on the determination of the upper endpoints of the spectra. 


These considerations lead to the following selection criteria, valid for all slepton studies: The events contain 
less than 10 charged particles, and two lepton candidates. The total charge should vanish, while each of the 
lepton-candidates should have opposite charge. The mass of each of the candidates should be less than M T . 
The total visible energy in the event, U vis , should not exceed 300 GeV, and the missing mass, M miss , should 
be larger than 200 GeV. No particle in the event should have momentum above 180 GeV. Lepton-candidates 
are found using the DELPHI r-finder 


18 , which is designed to identify both isolated electrons or muons, 


and decays of r leptons. It is also robust against extra activity in the detector coming from beam-beam 
effects and overlaid low-px 77 events. 


4.2.1 Selectrons and smuons 

As can be seen in Fig. [2] the selectron pair production cross section is huge in our scenario due to the i-channel 
neutralino exchange, allowing a very precise determination of the masses and polarised cross sections in a 
short time. At larger integrated luminosities, very precise measurements of the selectron and LSP masses can 
be obtained from the edges of the spectrum. In addition to the generic slepton selection criteria presented 
above, it is also demanded that both lepton candidates are identified as electrons. To separate (?r and eL, 
the beam-polarisation is helpful: for U+80,-30 > 6r production is enhanced, while eL production is suppressed. 
For 'P_8 o,+ 30) the situation is the opposite. In addition, since eL is appreciably heavier than than Hr, the 
decay products of the former are less boosted, and hence the detected leptons tend to be less back-to-back 
in eL events. A cut on the transverse momentum of one electron with respect to the other one, Pt 1 ^) is 
therefore applied to further separate the two cases: pti®* +Pr 2 el should be less than 30 GeV for eR candidates, 
and larger then 100 GeV for eL candidates. With these cuts, the selection efficiency is 51 % for eR and 47% 
for e L . 

Figure [9] shows the spectra of electron energies in selected di-electron events after collecting 500 ftV 1 of data 
for each of the beam-polarisations U_80,+30 and "P+80,-30- The Hr signal stands out above quite small SM 
and SUSY background for beam-polarisation "P+80,-30- The eL signal is less clean both due to smaller signal 
cross section (cf. Fig. [2]), larger background and lower efficiency, but is nevertheless quite prominent. In 
particular, both edges are detectable. This will still be the case even for substantially smaller branching 
ratios for the direct decay to lepton and LSP, compatible with LHC limits. 

The position of the edges is determined by sub-dividing the full data sample in sub-sets, and finding the 
most and least energetic lepton in each sub-set, after excluding a certain fraction of the extreme cases. The 
size of the sub-samples, and the fraction to be excluded is optimised to yield the lowest possible uncertainty 
on the endpoints. The resulting endpoints averaged over the sub-samples show a bias. This bias has been 
corrected for by means of a toy Monte Carlo procedure; the uncertainty on the SUSY masses is determined 
with this procedure as well. In this way, we obtain A-f^o = 95.47 ± 0.16 GeV and M-g r = 126.20 ± 0.21 GeV 
from the eR spectrum. The true masses in STCx are Af^o = 95.59 GeV and M-g R = 126.24 GeV. 

The large SUSY background to eL is mostly eReL productioi^] This channel gives important information 
on the neutralino mixing, since e.g. in case of light higgsinos the i-channel would be strongly suppressed by 
the small electron Yukawa coupling. In particular, if both beams are given left-handed polarisations, only 

3 Equation[2]shows that the lighter state receives a lower fraction of the total initial energy than the heavier state in eReL 
production, so that the upper edge from the eR decay is lower than what it is in the symmetric case in eReR production seen 
in Fig. [ 9 J 1 . 
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Figure 9: Electron spectra from selectron decays and background, with electrons from eR decays (a) and 
electrons from eL decays (b). 


the e + e —>• e R eL process is possible. As this reaction proceeds exclusively via neutralino exchange in the 
f-channel, its size gives insight to the neutralino mixing 78 . 
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Figure 10: Muon spectra from smuon decays and background. Fig. (a): muons from /tr decays and (b) 
muons from Ji l decays. 
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Figure [To] shows the spectra of muon energies in di-muon events, under the same conditions (/ Cdt = 500 ftU 1 
at each of the beam-polarisations). In this analysis, the generic selection is supplemented by demanding that 
both lepton candidates are identified as muons. The same criteria as for e are used to separate /Ir and ptL 
candidates. The selection efficiency is larger for g than for e: it is 65 % and 60 % for /1 r and /7 l, respectively. 
This is due to the t-channel contribution to e: even though the sleptons aren’t highly boosted, their decay 
products nevertheless carry the imprint of the initial slepton angular distribution. Therefore, the cuts 
designed to remove backgrounds at low angles to the beam-axis or with low missing transverse momentum 
tend to remove more e events than Ji ones. Using the same procedure to extract the edges, the LSP and 
smuon masses can be determined to M^o = 95.47 ± 0.38 GeV and M^ R = 126.10 ± 0.51 GeV, once again 
in good agreement with the true masses in the STCx model: M^o = 95.59 GeV and M^ R = 126.16 GeV. 
It can be noted that, as expected, the SUSY background to /IlMl production is much lower than for eL<3L 
production, and that a significant signal would be expected even if the branching ratio to the direct decay 
would only be a few percent. Combining the measurement of the LSP mass from the right-handed selectron 
and smuon analyses yields an uncertainty of = 147 MeV, ie. slightly above 1 permille. 
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Figure 11: Threshold scans at the e + e —► eReR (a) and e + e —>• /xr^r (b) thresholds. 


In addition to the mass determination from the spectrum edges, the mass of both eR and )Ir can be 
determined by scanning the production threshold near 250 GeV, as illustrated by Fig. 11 Close to threshold, 
the cross section is obviously small, but on the other hand, the signal is very clean: since the sleptons are 
produced almost at rest, and they undergo two-body decays, the decay products are almost mono-energetic. 
In addition, in STCx the mass-difference between the LSP and <3 r and )1 r is rather large, so that even 
in a decay at rest the produced leptons have momentum ~ 25 GeV. Hence, by selecting events with two 
opposite-sign, same-flavour leptons with large acoplanarity, and with momentum in a narrow window, a large 
efficiency, low background sample can be obtained, and an significant excess of events can be obtained quite 
close to the threshold. For STCx, it was demanded that the two leptons should have momentum between 
20 and 37 GeV, and both the acoplanarity and acollinearity angles should be below 3.1 radians. With these 
cuts, the efficiency for the signal is between 85 and 95%, and the signal-to-background ratio is above l for 
almost all points, the exception being the lowest £ cms for the / 7 r. 


Sleptons are scalars - if not the new physics the observations have revealed is not supersymmetry. Therefore, 
it is certain that slepton-pairs are produced in a P-wave, and hence that the rise of the cross section with 
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increasing E cms is proportional to /3 s = 


1 - 4 ( Mj/E c 


1 3/2 


. Investing a few months of ILC beam-time 4 

220 MeV, by fitting 


the mass of eR can be determined to ~ 190 MeV, and at the same time that of /7 r to 
the observed background-subtracted cross section to (3 3 (E cms ). 

In addition, one can test the hypothesis that the observed states are spin-1/2 particles, rather than scalars. 
If the particles indeed have spin-1/2, the pairs are produced in an S-wave, and the rise of the cross section 
with increasing E cms would be proportional to /3 rather than /3 3 . The dashed curve in Fig. 11 is the best fit 
of /3(E C ms) to the data. For the Hr case the fit-probability is < 10 -9 , while it is 7 x 10 -5 in the /Ir case. 
The spin-1/2 hypothesis would therefore be excluded by the data. 


4.2.2 The r-sector 


Especially in r-coannihilation scenarios, a precise determination of the r sector is essential in order to be 
able to predict the expected relic density with sufficient precision to test whether the \i is indeed the 
dominant Dark Matter constituent. The ri-pair production is different from eR- or /tR-pair production in 
several aspects: The mass difference to the LSP is much smaller, meaning that the r spectrum is softer 
than the spectrum of the leptons in eR or /tR decays. In addition, r leptons decay, further softening 
the spectrum of observed particles, and making the particle identification requirements less effective in 
background suppression. This leads to a signal that much more resembles that of 77 events, but also more 
resembles di-boson events decaying to tvtv. The generic slepton selection therefore needs to be supplemented 
by several further criteria to reduce these sources of background: The requirements on E v ; s and M m ; ss are 
strengthened to < 120 and > 250 GeV, respectively, and the visible mass, M vis , should be below M z — 5, 
which reduces the di-boson background. The cosine of the direction of the missing momentum is required 
to be between -0.8 and 0.8, M v ; s should be above 20 GeV, and the total energy observed below 30 degrees 
to the beam-axis should not exceed 2 GeV. This selection reduces the 77 background, which is then further 
decreased by a cut on the likelihood that the event is a 77 event. Finally, to reduce the SUSY background 
from eR- or pR-pair production, as well as from di-boson events, it is required that the event is not identified 
as a di-electron or di-muon event. With these cuts, the selection efficiency for Ty-pair production is 17%. 


Only the upper endpoint can be measured in t\ production: due to the decay of the r, the lower endpoint is 
only visible as a knee in the spectrum of the decay-products of the r. Because of the small mass-difference 
between the T\ and Xi, this knee is in a region where the spectrum is strongly distorted by the cuts removing 
the 77 background. Contrary to the case for e or /t, the upper kinematic limit is not an edge, but the 
endpoint of the spectrum of r decay-products. This endpoint is determined by fitting the background in 
the region well above the endpoint and then fitting a signal contribution in the data above the extrapolated 
background fit. 


Figure [T2] shows the energy-spectrum at E cms = 500 GeV of selected r-jets for an integrated luminosity of 
500fb _T 7and polarisation V+so,- 30 - In the case of 77 production, the endpoint could be determined to be 
Eendpoint = GeV, corresponding to an uncertainty of of 200 MeV, if an uncertainty on the LSP 

mass of ~ 100 MeV is assumed. Similarly, the t/ mass could be determined with an uncertainty of 5 GeV. 


In 75 , where a model quite similar to STCx has been studied, it is found that, in addition to , the 


production cross section for both these modes can be determined at the level of 4%, and the polarisation of 
r-leptons from the r-\ decay, which gives access to the t and \i mixing could be measured with an accuracy 
better than 10 % from shape of the 7r spectrum in the r -7 7 t + v t mode or to better than 5 % by a template 
fit of R = En/Ej et in the t~ — > p~v T — > (and c.c.) mode. 


4 Note that this centre-of-mass energy range is also optimal for studying the properties of the higgs boson with the model- 
independent recoil-mass method. 

'"'Interaction of sfermions and gauginos conserve chirality, while the Yukawa interaction of the higgsinos flips chirality. 
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Figure 12: The r spectra in t\ decays and SM as well as SUSY background, with an endpoint fit (a), and 
t spectra in r 2 decays and background, with endpoint fit (b). 


4.2.3 The sneutrinos 


While the vast majority of sneutrino decays proceeds to a completely invisible final state, the v T has a 
branching fraction of about 5% to T\W. In this situation there are two possible strategies: a) search for the 
completely invisible final state via its recoil against a hard photon from initial state radiation in analogy to 
Section 4.1.3 and b) select v T pair events where at least one visible decay occurs. Neither case has yet been 


studied in detailed simulations, but we will sketch the strategies here. 

In case of the ISR recoil, the situation is more difficult than for xi pair production since the cross section for 
the low-background polarisation U+80,-30 is about one order of magnitude smaller than for the xi case. At 
the same time, due to the larger mass of the sneutrinos, the energies of the photons are smaller, thus buried 
in the steep shoulder of the x? spectrum 
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Therefore, the 'P_8o,-t-30 combination seems more promising. 
In this case, the sneutrino pair production cross section is an order of magnitude larger than the cross section 
of Xi pair production, and of roughly the same size as the xi cross section in the other polarisation, c.f. 
Table [7] The price to pay is a background from SM neutrino pair production which is about a factor of 5 
larger than in the 'P-1-80,-30 case 76 . This means that the control of systematic uncertainties becomes even 


more important than in the classic xi case. However, due to the lower photon energy endpoint, a large part 
of the photon energy spectrum will be signal free, which should give an excellent possibility to constrain 
absolute normalisation uncertainties as well as shape uncertainties, arising e.g. due to the finite knowledge 
of the beam energy spectrum. Therefore, based on the experience from radiative neutralino production, we 
expect that both cross-section and mass measurements are possible, but would need a quantitative study. 

The visible decays of the v T lead to a rather unique final state with a r, a (hadronic) W and large missing 
four-momentum. SM background from WW —> rv T qq' as well as the main SUSY background from xfxf —> 
r-| v T xiW can be reduced very effectively by choosing U+8 o,- 30 i since the polarisation dependence for the 
pure s-channel v T is much weaker. Furthermore, both backgrounds feature different kinematics from the 
signal, since the r and the W stem from different hemispheres, and not from the same as in the signal 
case, which can be exploited efficiently once the masses of the fy and xf have been measured. The SM 
WW background can be further suppressed by exploiting the forward-backward asymmetry 75 . Once a 
sufficiently clean signal has been selected, the v T mass can be determined from the endpoints of the W 
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energy spectrum (c.f. Section 4.1.2), based on the ri mass measured in direct ti pair production as described 
above. Similar to the ISR recoil case, this analysis seems feasible, but awaits a detailed simulation study for 
quantitative conclusions. 


Thus, the currently most obvious way to access the v mass are the Xi cascade decays discussed in Sec¬ 
tion [TL2J 


4.3 Analysis of sleptons in cascade-decays 


A particularly interesting channel for the determination of slepton properties is e T e 


+ o” 


X2X2 and the X° 

decay to mr. g (or equivalently to eRe), even if the branching ratio is at the level of a few percent as in 
STCx. These cascade decays can be fully kinematically constrained at the ILC, and would promise to yield 
even lower uncertainties on the mr and Hr masses than the threshold scans, of the order of 25 MeV. This 
is estimated on an earlier study in a scenario with about twice as large branching ratios for the considered 
decay mode, where a precision of 10 MeV 80 was found. The corresponding distribution of the reconstructed 
/Ir mass is shown in Fig. [Id}' a), including all SM and SUSY backgrounds. 


Also for the study of the ti, the channel e + e - —tX2.X2 and the X 2 decay to t\t is quite powerful. This 
decay channel has a much larger branching-ratio than the electron and muon channels, but the presence 
of the neutrinos in the decay of the r prevent the exact kinematic reconstruction of the events. However, 
approximate reconstruction is still possible, as shown in Fig. |13| (b). Potentially, this method could yield 
comparable results to a threshold scan. 




Figure 13: Reconstruction of slepton masses from X 2 X 2 — ► Mil i n SPSla, which has a very similar spectrum 
to our case. We show the reconstructed Mr mass (a) and the reconstructed T\ mass (b). From 
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5 LHC-ILC Interplay 


In this section, we employ the simulation studies based on the STCx scenarios to illustrate how discoveries 
and measurements at the LHC and a future linear collider like the ILC could work together to gain as precise 
knowledge as possible about the origin of the beyond-the-Standard-Model (BSM) observations. 
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5.1 Discoveries at the LHC 


By construction, several sectors of the STCx spectra offer discovery opportunities at the LHC. For instance 
the left-handed selectron and srnuon could be seen early if their branching ratio for the direct decay tj, —»• Ixi 
is large, c.f. Section |2.1[ The right-handed sleptons, on the other hand, are much harder to detect at the 
LHC, due to the - for the same mass values - lower production rates, and in particular in the STCx points 
due to their smaller mass difference, which leads to softer leptons in the visible final state. The t\ and T 2 
will be very hard to observe. 


Figure [14] summarises the significances obtained by the analyses presented in Section [3] as a function of the 
integrated luminosity and as a function of the assumed systematic uncertainty on the background estimate. 
The first significant deviation of the coloured sector of the STCx scenarios from the SM would appear in 
the single-lepton channel at the LHC, caused by a mix of top- and bottom-squark pair production, as well 
as squark-gluino production. Assuming a systematic uncertainty of 15 % on the background prediction, a 
5 a deviation could be visible with less than 50 fb” 1 . The amount of integrated luminosity required rises 
strongly with the masses of the top and bottom squarks, as can be seen by comparing Fig. [l4|a) with the 


lower top and bottom-squark masses in model STC8 to Fig. 14 
STC10 with 200 GeV higher top- 


c), where the same is shown for the model 


and bottom-squark masses. From the existing LHC analyses, one would 
not necessarily expect a specialised search for direct top-squark production to be one of the first discovery 
channels. However, the limits on direct top-squark production are usually determined for one or maybe 
two decay channels, assuming that direct top-squark pair production is the only production mechanism of 
SUSY particles, while in the full models there are several production and decay channels contributing to the 
discovered signal. In our case, we see also decays of top and bottom squarks to final states with higher-mass 
neutralinos and charginos, and contributions from gluino-pair production. If, however, the gluino mass would 
be higher, the direct top- and bottom-squark production would be the main contribution, leading to a lower 
significance and later discovery. 


In case of STC8, the dedicated search for bottom squarks decaying directly into a b-quark jet and the 
LSP would report a 5er deviation from the SM with about 90 fb _1 (assuming a systematic uncertainty of 
15%). With the higher bottom-squark mass of STC10, however, a substantially larger amount of integrated 
luminosity, 300fb _ , would be required. This is due to the fact that the bottom-squark selection has a much 
higher purity, in particular it is less sensitive to gluino-squark production than the top-squark search. Thus 
the difference in integrated luminosity required to observe a signal in these channels does not necessarily 
correlate with the masses of the targeted squarks. 


After an integrated luminosity of 150 fb _ , the multi-lepton channel would report a 5 a deviation from the SM 
in either STC8 or STC10. Since b-tagging does not play a role in this analysis, it will be highly likely that this 
deviation is not caused by third-generation partners, but by lighter sparticles, i.e. electroweakinos and/or 
sleptons. The variation of the excess over the numerous signal regions gives a hint that several different 
sparticles with sizeable mass differences contribute to the overall signal. For the expected significance shown 
in Fig. 14 we combine the results of all search regions 81-83 . The prediction for the non-pronrpt background 
is the most controversial part of this analysis, since the fake rate has to be assumed. With the fake rate 
increased by a factor of 2, the discovery would require an integrated luminosity of 200 ftU 1 instead of 
150 ftT 1 . 


A further hint towards the existence of rather light third-generation partners would come a little later 
from the inclusive hadronic search: While this channel would remain at the 2 to 3 cr-level before the High- 
Luminosity LHC when no b-tagging is applied, this could be enhanced in case of STC8 to the 5 er-level after 
less than 200 fbC 1 by requiring two jets to be identified as b-quark jets, as shown in Fig 14 In STC10, 


about 600 fb~ would be needed for a discovery in this channel. In addition, it will be extremely difficult 
to separate the contributions from direct top- and bottom-squark production and gluino decays via the 
third-generation squarks in this analysis. 


All search channels show strong dependency on systematical uncertainties at different pile-up and higher 
luminosity scenarios, as shown in Fig.[l4](b) and (d). In particular the single-lepton stop (Section[T3|) and the 
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Figure 14: LHC searches: Expected significance as a function of luminosity and of the systematic uncer¬ 
tainty, respectively, for the two signal scenarios STC8 (a,b) and STC10 (c,d). For plot (a) and (c) we assume 
a systematic uncertainty of 15 % for the bottom- and top-squark search, and 20 % for the other searches. 
As the spectrum of the electroweak sparticles is the same for STC8 and STC10, we show the multilepton 
analysis which focuses on the electroweakinos only for STC8. 


hadronic searches would profit significantly from an improved understanding of the background level beyond 
the 15 — 20 % assumed as default here. The hadronic search requiring two 6-tagged jets (Section [iO] ) and the 
hadronic bottom-squark search (Section 3.1l are most susceptible to systematic uncertainties. For STC8- 
like masses, they would miss the 5 a -level with 300 fb -1 if the systematic uncertainty on the background 
was larger than 25% and 30%, respectively, while in case of STC10 a 5er discovery is only possible if the 
backgrounds are controlled to better than 15 %. The single-lepton search is, not surprisingly, more robust and 
would still reach 5 cr with uncertainties of 50 % in case of STC8 and 40 % in case of STC10. The sensitivity of 
the multilcpton search depends on the significance the 45 different search bins. With decreasing systematic 
uncertainty, the contribution of the low mr signal regions gains significance, leading a steeper rise in the 
sensitivity can be seen for very small systematic uncertainties, as shown in Fig. |14[ c). However, we expect 
systematic uncertainties of the order of 20 % for this analysis. 
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5.2 Signal characterisation at the LHC 


Once a clear deviation from the SM has been discovered, the immediate question will be the origin of this 
deviation. Given the fact that at this stage neither the masses of the produced sparticles, nor their decay 
chains and branching ratios are known, it seems highly unlikely that the contributing production modes 
could be identified at that stage. 

The most promising candidates for isolation of a sufficiently pure sample of an individual decay chain are 
the direct slepton decays l\\ and the direct bottom-squark decay bi —> fry?- These could then be used 
to obtain information on the masses of the produced sparticles via kinematic edges. 

For instance one can use the so-called boost-corrected contransverse mass (toct) [29j| 30 , which is defined as 

to ctC?i> 32 ) = [£t0i) + E T {j 2 )] 2 - [pt(ji) - PtUi)] 2 

= 2 p T (ji)PT(j 2 )(l +cosA0(ji, j 2 )) , 


can be used for extracting information from events in which two heavy particles decay into a visible part jj 
and missing energy. This mass is invariant under equal and opposite boosts of the parent particles in the 
transverse plane. For parent particles produced with small transverse boosts, mcT is bound from above by 
an analytic combination of particle masses. 

In the case of the direct bottom-squark decay, the distribution of the contransverse mass of the two jets, 
toct (ji; J 2 )! is expected to show an edge at 


*c# e = ^ m l°)A 


(4) 


which equals 780 GeV for STC8, and 1000 GeV for STC10. 


The expected distributions are shown in Fig. 15 for an integrated luminosity of 300 fb after the selection 
presented in Section [372} With this amount of integrated luminosity, the edge position in model STC8 could 
be determined, while for STC10 this would be possible only with the luminosity that can be reached at the 
HL-LHC, 3000 fb 


-1 


The edge is smeared out by the jet energy resolution, which has still to be determined for the 14TeV data 
with the corresponding jet reconstruction at higher pileup. Assuming a resolution of 10% for PtU) = 
300 — 500 GeV jets (which is a typical transverse momentum for jets from bottom-squark decays, as shown in 
Fig-0- motivated by a corresponding 7TeV measurement [84], we can deduce a resolution of about 15 % for 
the Toct variable, corresponding to about 120 GeV in STC8 and 150 GeV in STC10. The edge position gets 
further diluted or even biased by the SM background as well as by the SUSY background, in combination 
with fluctuations due to low statistics. While the situation could be partly improved by removing the SM 
background based on prediction from simulation, which is expected to be well understood with 300 fb -1 of 
data, the SUSY background remains as a priori unknown and can distort the determined endpoint as shown 
in the following. 

In order to determine the mass edge, one could e.g. exploit a typical edge-finder like the so-called edge- 
to-bump method 85 . To test the capability of this method for the model STC8, we generate 1000 times 


pseudo-data from the distribution shown in Fig. [T5[a), and determine the edge position with the edge-to- 
bump method. The mean of the fitted value is at 832 GeV with an RMS of 114 GeV, which is compatible 
with the expected value of 780 GeV. This RMS is expected from the resolution of the toct variable. With 
ten times more statistics from the HL-LHC, this uncertainty could be reduced roughly by a factor of three. 
Without further knowledge about the Xi> the determined endpoint can be interpreted as lower limit on the 
bottom-squark mass using Equation [4j assuming a mass-less neutralino, of 832 ± 114 GeV. Higher masses of 
the neutralino would lead to higher bottom-squark masses. 

Analogously, an edge in the contransverse mass spectrum would be expected from the decays of eb and /!l 
near rn^ e (Zi, I 2 ) = 175GeV. Without any information from ILC, this would yield a lower limit on the 
slepton mass. 
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Figure 15: LHC bottom-squark search: Comparison of the mcT distribution for signal and background 
from SM and other SUSY processes for STC8 (a) and STC10 (b). All SUSY processes from the STCx models 
are considered, with the signal being stacked on top of the stacked background. 


5.3 Discoveries at the ILC 

The first channel to manifest itself at the ILC depends on the assumed running scenario. If the ILC starts 
out as a higgs factory at E cms = 250 GeV, then e + e“ —> T\T\ and XiXil would be the first observable 
channels, while e/j- and pair production is just beyond reach. The measurement of the r mass, however, 
would be challenging close to threshold, since both upper and lower edge of the r-lepton energy spectrum 
would be in the region affected by background from multi-peripheral two-photon processes. 



EjJGeV] 


Figure 16: Momentum spectrum of events with e + e— and missing 4-momentum. The assumed luminosity 
of Sfb^ 1 corresponds to one week data taking at design luminosity. 

On the other hand, the LSP mass and pair-production cross section could be measured at least with a 
few percent precision from the energy (or recoil mass) spectrum of the accompanying initial state radiation 
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photons as discussed in Section |4.1.3| Although that study has been performed at a higher centre-of-mass 
energy, the quality of the mass determination is expected to be comparable at E cms = 250 GeV, since the 
photon energy spectrum less spread out while the size of the cross section is similar. Since the neutralino 
pair-production is dominated by t-channel selectron exchange, the mass of the lighter selectron and its 
helicity can be determined from the measurement of the polarised cross sections. However, this requires a 
sizeable amount of integrated luminosity collected with right-handed electrons and left-handed positrons in 
order to control the SM neutrino background. 


In the recently published running scenario 23 , the ILC starts up in its TDR baseline configuration at 
E C ms = 500 GeV. In this case, the picture changes drastically with already a modest amount of data, 
since very first evidence for BSM at the ILC would be observed after a few days of running from selectron 
pair production. Figure 16 shows the selectron signal expected with only 5fb _ , allowing already a first 
determination of both the selectron and LSP masses according to the technique described in Section [4] 
Since the lighter set of sleptons must have a 100 % branching ratio to lepton and LSP, their production 
cross sections can be measured unambiguously. In particular for the selectron and smuon case, where the 
mixing is typically small, the coupling can be extracted once the mass is known, thus enabling the first 
verification that the couplings of the exotic states equal those of their SM partners - a fundamental property 
of supersymmetry! 


Beyond the specific case of STCx, any sparticle lighter than half the centre-of-mass energy of an e + e“ 
collider will be observable. If furthermore, like in the STCx models, i?-parity is conserved so that the LSP 
is stable, the NLSP will have a 100% BR to the LSP and the corresponding SM particle. This means that 
the search for NLSP pair production at the ILC can be seen as a loop-hole free search for SUSY 86 . The 


case of a general STC-like model, where the NLSP is the ti, is illustrated in Fig. 17 showing the exclusion 
and discovery reach for the NLSP at the ILC is shown in the plane of the only free parameters, i.e. o and 
Mnlsp • It should be pointed out that the case with the NLSP being the t\ is one of the most experimentally 
difficult ones, due to the need to detect r-leptons in the final state. 




Figure 1 7: Discovery-reach for a t\ NLSP after collecting 500 fb 1 at U cms 
full scale, Right: zoom to last few GeV before the kinematic limit. From 86 


500 GeV at the ILC. Left: 
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5.4 Signal characterisation at the ILC 


23 


500 fb would be collected at 


According to the currently favoured running scenario for the ILC 
E C ms = 500 GeV within the first 4 years of operation. For these data, a sharing between the four possible 
configurations with different polarisation signs is foreseen: 40 % of the data would be collected at each of the 
configurations V-so ,+30 and V+so,- 30 , while 10 % would be collected at each of the V-so ,-30 and "P+ 80,+30 
ones. This would allow a first assessment of all SUSY particles with masses below 250 GeV, including their 
masses and mixings, with a precision typically a factor 1.6 worse than then values quoted in section [4] 

After the initial run at E cms = 500 GeV, it is foreseen the lower the centre-of-mass energy to scan the top-pair 
production threshold and to run near the Zh threshold for a high-precision determination of the higgs boson 
mass and its coupling to the Z boson. These runs would also be of high interest for SUSY spectroscopy: 

Running near the Zh threshold could include a scan of the threshold for pair production of the lighter sleptons 
as shown in Fig. EH which provides mass measurements of the en and with precisions of 190 MeV and 
220 MeV, respectively. The shape of the threshold as well as angular distributions identify the produced 


partners of electrons and muon as scalar particles 87 


At or slightly above the tt threshold, the cross section for T 1 T 2 mixed production could be measured with 
high precision due to the absence of the r 2 -pair production background which limits the purity of this 
measurement at E cins = 500 GeV. This provides an interesting possibility to determine the t mixing and 
tan/3. Together with the masses of the T\ and the and the Xi mixing determined at E cms = 500 GeV 


(c.f. Section 4.2.2 and 4.2.1), these are important inputs for the prediction of the dark matter relic density 


within the MSSM, since in the STCx scenarios LSP pair annihilation and ri-coannihilation contribute about 
equally to the cosmic dark matter annihilation. By comparing the predicted value for the relic density to 
the cosmologically observed one, the Xi can be either identified as the sole constituent of dark matter, or as 
being responsible for only a fraction of the observed relic density. This has been explicitly demonstrated in 
a SUSY scenario with an electroweak sector very similar to the STCx cases 


After a luminosity upgrade, the integrated luminosity at 500 GeV is then foreseen to be increased to 4ab , 
with the same polarisation sharing as before, enabling a full BSM precision program, which reaches the 
permille-level for many observables. While the currently assumed running scenario is based on guaranteed 
measurements only, it will be adjusted to future discoveries by scheduling further threshold scans as well 
as runs at dedicated energies not far above important thresholds. Such special runs give e.g. important 
information for parameter determination in the neutralino sector, which in turn enables predictions of the 
masses of the kinematically not accessible states 89 . The results of the initial run at E cms = 500 GeV will 


give decisive input to determining the detailed running strategy. 


5.5 Combining ILC and LHC 


With ILC information on the lower part of the spectrum, the situation at the LHC changes drastically. While 
the clear advantage of the LHC is the larger kinematic reach for the production of the more heavy sparticles, 
disentangling their signals is challenging due to the multitude of decay modes, most of them with small 
branching ratios. Thus the information from the ILC might provide the key to access all the information 
contained in the LHC data, on one hand by exploiting the knowledge of the lower parts of the decay chain 
in the data analysis, on the other by combining LHC and ILC results in global pMSSM fits. This has been 
pointed out in several previous studies 88,90,91 . Most of them, however, are based on scenarios with a 


light, sub-TeV, coloured sector, which is by now excluded by LHC measurements. Thus we highlight here 
the specific possibilities arising in the more challenging STCx benchmark points. 


First of all, the determined edge position from bottom-squark pair production could be turned into a mea¬ 
surement of the bottom-squark mass by adding the information from ILC about the Xi mass, which can 
be determined at the ILC with a negligible uncertainty. Using Equation [4j we find a bottom-squark mass 
of 843 ± 115 GeV, which is compatible with the true value of 795 GeV, as well as with the value obtained 
for a mass-less LSP. Even with HL-LHC precision of 30 — 40 GeV, the hypothesis of a mass-less LSP would 
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yield a bottom-squark mass compatible with the model value. However, already for an LSP mass of 200 GeV 
which could still be accessed at the ILC with E cms = 500 GeV from radiative neutralino production (c.f. 
Section 4.1.3), the bias due to the zero-mass hypothesis rises to 50 GeV and thus surpasses the statistical 
uncertainty on rn^ e . 


In addition, we could use the measured bottom-squark mass to calculate the cross section for this process, 
and compare it to the cross section times branching fraction determined above from counting the events in 
the signal region. While the mass determination is independent of the spin of the newly discovered particle, 
the measured event rate that can be related to the cross section (times branching fraction) depends on 
the spin of the new particle and therefore help identify whether the newly discovered particle is actually a 
SUSY particle or not. Final conclusions could only be drawn if the branching fraction can be determined 
in addition. For the latter, input from the ILC on the lower-mass sparticle spectrum might be exploited as 
well. 


In the electroweakino sector, since now not only the masses, but also the mixings and decay modes of the 
lighter electroweakinos are known, they enable a significant improvement of the analyses of LHC data. In 
case of the multi-lepton analysis, the signal contribution from the lighter electroweakinos can be subtracted, 
giving access to the heavier electroweakinos not accessible at the ILC with E cms = 500 GeV. 

One example is given in the following. We want to study the decay of a > Xi Z —> v £ + £~. Therefore, 
we select events with three leptons and a large amount of missing transverse energy. The opposite-sign same 
flavour lepton pair £ + £~ covers the full information about the Z-boson and can be identified by requiring the 
invariant mass as close as possible to the Z-boson mass. In the following we refer to the remaining lepton 
from the xt decay when mentioning the lepton, and assume that the masses of the Xi and v will be known 
with sufficient precision from ILC measurements, as discussed in section |4.1.2| The invisible decay products 
of the v lead to signatures with sizeable missing transverse energy. Unfortunately this is not the only source 
of missing transverse energy, since the decay of the other initial sparticles will also produce at least one LSP, 
hiding the information about the momentum of the v. In order to avoid this, we boost into the Xi rest 
frame, where the lepton has a fixed energy E’ which can be calculated by requiring energy conservation: 


E 'lep ~ 


( A %) 2 


+ ( miep ) 2 - ( top ) 2 


2M xf 


8.6 GeV 


(5) 


with the masses known from ILC: M~± =206.1 GeV, ?np=197.3 GeV and mi ep « OGeV. The transformation 
from the lab frame into the xt rest frame is done via the boost vector f3: 

E lep = l E lep - 7/3 • Plep, ( 6 ) 


with 7 = 1/y 1 — /3 2 . This vector has three degrees of freedom, the magnitude and two angles. For the 
angles we assume that the boost is parallel to the lepton. This is a good approximation in the most cases, 
since the Xi decay does not add much to the boost of the lepton. With this assumption we can calculate 7 : 


1 / E lep Elep\ 

2 \E lep E[ ep ) 


(7) 


and have the full information about the boost vector (3. As next step we boost the Xi back into the lab 
frame (P~± = +7/3 M~± ). As last step we add the four-momenta of the Z boson and the xt in order to 

Xi Xi 

calculate the invariant mass of x 2 • 


Figure 


18 


shows the the M~± by using the MC level information as well as reconstructed objects. The 


resolution gets smeared by two different effects: First, the E[ is not fixed because of the natural width of 
Second, the boost vector /3 is not perfect parallel to the lepton. Therefore, we expect to 


the M~± and to 

Xi 


measure the M~± to be 412 ± 43 GeV. In order to suppress the SM background we require exactly three 
leptons (jjt = 25/15/10GeV), one opposite sign same flavour lepton pair with an invariant mass between 
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84 GeV and 96 GeV, at least one b-tagged jet (pt > 100 GeV and p < 2.4), more than three jets (p r > 40 GeV 
and p < 2.4) where the leading jet has px > 120 GeV and U™ ISS > 200 GeV. With this selection our signal 
consists dominantly of \f produced by stop decays, which has the best signal to background ratio. With 
a slightly softer selection we would also have sensitivity to direct xf Xm production, which also can be 
accessible at HL-LHC. Per construction the M~± must be larger than ./(Af~±) 2 + ( Mz) 2 > 200 GeV. As 

first observation we see much more SUSY events than we would expect from our targeted decay chain. Most 
of these events include heavy electroweakinos (xf , x®, x°) decaying via Z bosons. The third lepton steins 
from the other decay products, which in most cases are r leptons. This signature carry mass information 
about the heavy charginos and would be worth to be studied with the hadronic r lepton final states. The 
other SUSY events passing this selection contains slepton decays, which has a kinematic edge signature close 
to the Z mass. This kind of background could be studied further if one requires the invariant mass of the 
opposite-sing same-flavour lepton pair exactly at the position of the edge. Overall the peak from our targeted 
decay chain around 420 GeV is visible. 
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Figure 18: xf m ass information at LHC. The left plot (a) shows the MC truth information of the invariant 
mass obtained with the procedure described in the text as well the reconstructed M~± for those events. The 

A2 

right plot (b) shows the reconstructed M~± for all events we would expect. There the SM and the SUSY 
events are stacked together. 


The obtained knowledge of the xf an d X 3 4 masses in turn would add significantly to the physics case of a 1- 
TeV-upgrade of the ILC, or of an even higher energy e + e _ collider like CLIC. Furthermore, the full knowledge 
of the electroweakino masses, mixings and decay modes could provide the decisive information to isolate a 
signal from top-squark production, e.g. in it’s largest decay mode U —> b\f, enabling a determination of the 
ti mass. 


6 Conclusions 


In this paper we discussed the complementarity and interplay of a proton-proton collider, the LHC, and 
and electron-positron collider, the ILC, in discovering new particles and in determining their properties. As 
example we used an i?-parity and CP conserving supersymmetric model where the f) is the next-to-liglitest 
supersymmetric particle and has a small mass difference of about 10 GeV to the lightest supersymmetric 
particle, the LSP, which is the lightest neutralino \i- I n such a scenario, ri-coannihilation allows for a suffi¬ 
ciently small dark matter relic density. All sleptons and electroweakinos have masses below 500 GeV, while 
the lightest coloured sparticles, the lighter top and bottom squarks, have masses around 800 or 1000 GeV. 
All other coloured sparticles are much heavier, up to 2 TeV. 
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We showed that such a scenario can be easily discovered at the LHC running at 13/14 TeV. In particular, 
this is true for the heavier selectrons and snruons, the heavier electroweakinos, as well as the lighter top and 
bottom squarks, if their masses are not much higher than a TeV. Depending on the search channel, on the 
exact masses and on the achieved control of systematic uncertainties, the integrated luminosities to discover 
deviations from the SM expectation due to production of these sparticles range from 50 and 1000 fb -1 . The 
earliest discovery would come from the single-lepton stop search, followed by hadronic searches including 
b-tags. Inclusive hadronic searches would require systematic uncertainties to be controlled better than 10% 
in order to achieve a 5 a discovery. 

However, in most cases the observed deviations from the SM expectation cannot be attributed to a single 
process, but result from a mixture of e.g. diverse electroweakino processes or, in case of the coloured sector, 
of stop, sbottom and gluino production. The best chances are to identify a single process with sufficiently 
high purity, in order to learn something about the properties of the produced sparticles originating from 
direct decays of the produced sparticle to its standard-model partner and the LSP. For the investigated 
scenarios, the direct decay has a sizeable branching fraction in case of the heavier sleptons and the lighter 
bottom squark. Here, a kinematic edge can be isolated e.g. in the contransverse mass distribution. Without 
further knowledge of the mass of the LSP, the position of this edge can be converted into a lower limit on 
the sparticle mass. 

As example, we analysed the case of bi —»• by^, and found that the edge in the contransverse mass distribution 
for the model with bottom squark mass could be determined with an uncertainty of about 115 GeV with an 
integrated luminosity of 300fb -1 at 14 TeV. Assuming that systematic uncertainties play only a minor role, 
this could be improved to about 30 GeV at the High-Luminosity LHC. For higher bottom squark masses 
of about 1 TeV, the discovery is still possible, however the uncertainties on the mass determination will be 
significantly larger due to strongly reduced number of signal events. 

The ILC on the other hand would complement these spectacular discoveries at the LHC with a systematic 
precision analysis of the lower part of the spectrum, actually discovering some of the lighter states. In this 
paper, we especially present up-to-date detector simulation studies of the slepton sector at the ILC, both in 
the continuum and in threshold scans. In particular a sufficiently light selectron will lead to a striking signal 
within a few weeks of ILC operation. Already after the 500 fb -1 collected in the first 4 years of ILC operation, 
this would allow a permille-level determination of the LSP mass, and the masses of the right-handed sleptons, 
improving further with the 4ab~ : foreseen for the full ILC program. In conjunction with threshold scans and 
operation with different beam polarisations, all sleptons and the lighter half of the electroweakino spectrum 
can be characterised with at least percent-level precision, including masses and mixing angles. This would 
demonstrate that the observed new particles are indeed of supersymmetric nature, and at least some of the 
parameters of the underlying SUSY model can be determined. 

This information from the ILC then creates new opportunities to analyse the LHC data. Obviously, a 
precise determination of the LSP mass from the ILC can be employed to turn edge determinations into 
mass measurements. Even better, the detailed information from the ILC can be used in the analysis of the 
data themselves to disentangle the contributions of different production modes and to reconstruct quantities 
which are sensitive to masses of the heavier sparticles with more complex decay chains. As an example, we 
illustrated this in case of the xt • With the knowledge of the masses of the \\ and the v from the ILC, 
the signal can be isolated from the electroweakino mix and its mass can be reconstructed on an event-by¬ 
event basis with a resolution of about 50 GeV with an integrated luminosity of 300 fb -1 at 14 TeV, with 
corresponding improvements at the High-Luminosity LHC. 

Although we studied the capabilities and the interplay of LHC and ILC based on a specific example, many 
aspects are transferable to other scenarios which comprise new particles in the kinematic reach both colliders. 
In this context, it should be noted that the LHC in many cases the sensitivity to the lightest new physics 
states is significantly smaller than to some of the more heavier states, and that some of the lighter states 
might even await explicit discovery at a lepton collider. We finally conclude that the combination of LHC 
and ILC data could reveal significantly more information about the properties and the origin of new particles 
than the results from either collider alone. 
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7 Appendix 


Outflow tables and numbers of events in the signal regions for the different LHC searches are given in this 
appendix. The outflow and event yields of the SM backgrounds and the signal samples for several signal 
regions in the hadronic search are shown in Table [9] 

Table 9: LHC all-hadronic inclusive search: Background and signal event yields corresponding to 300/fb. 
The vector boson denoted by ”V” refers to W, Z and 7 


Cutflow 

tt +jets 

V+jets 

VV+jets 

Top + jets 

Total SM 

STC8 

STC10 

Preselection 

3.93xlO Y 

3.94xl0 8 

4.5x10* 

5.7x10* 

4.43x10* 

27695 

11510 

HT > 1000 GeV 

1 . 6 xl 0 6 

1.04xl0 7 

236392 

120527 

12.3xl0 6 

8261 

3246 

MHT > 500 GeV 

8503 

108477 

4795 

483 

122458 

3502 

1895 

\A<j>(Jets,MHT)\ 

1094 

22256 

934 

43 

24328 

912 

519 


The outflow with the number of expected signal and background events for the bottom squark search is given 
in Table ITOl 

Table 10: LHC bottom squark search: Event yields for the different SM processes and STC models after 
applying each cut, the numbers are normalised to total 300fb _1 . 


Cutflow 

VV+jets 

tt +jets 

V+jets 

Top+jets 

Total SM 

STC8 

STC10 

Preselection 

86940500 

173553000 

11212700000 

34843100 

11508036600 

1160340 

1202280 

Lepton(e,^t)veto 

63661400 

86844100 

9884430000 

24245100 

10059180600 

936391 

974139 

N-Jets >=2 

19224300 

77540200 

3895480000 

16818900 

4009063400 

191599 

185023 

l st Jetp T >300 GeV 

688728 

2606200 

42888600 

388099 

46571627 

27418 

18690 

2” d Jetp T >200 GeV 

437937 

1602890 

31104200 

237408 

33382435 

15593 

9458 

2 nd Jetpr <70 GeV 

124572 

207703 

9029620 

60479 

9422374 

2801 

2130 

N-bjets =2 

543 

29836 

32045 

3290 

65715 

521 

159 

E™ 1SS > 450 GeV 

5.9 

386 

72 

14 

478 

214 

82 

to t > 500 GeV 

4.1 

12 

51 

0.27 

67 

184 

72 

H t > 850 GeV 

1.4 

2.8 

24 

0.05 

28 

74 

44 


Table [Tl] shows the expected number of SM background and SUSY signal events after each step of the event 
selection for the top squark search. 

Table 11: LHC top squark search: The event yields for the inclusive signal samples and several SM processes 
with 300 ftV 1 at 14TeV with 50 pileup interactions. 


Cutflow 

VV+jets 

Cutflow 

VV+jets s 

Top+jets 

Total SM 

STC8 

STC10 

Exactly 1 lepton 

10141300 

43887300 

548959000 

7665270 

610652870 

114405 

112180 

Njets >= 5 

57851 

3298680 

827110 

29067 

4212708 

4149 

1911 

1 or 2 b jets 

19249 

2405190 

232049 

21492 

2677980 

2642 

1040 

E^ 1SS > 400 GeV 

219 

5339 

1601 

62 

7222 

642 

410 

A cj) > 0.8 

194 

3990 

1397 

45 

5629 

536 

338 

Centrality >0.6 

117 

2275 

822 

24 

3240 

418 

282 

to t > 260 GeV 

4 

94 

0 

0 

98 

194 

134 

2 > 260 GeV 

2 

26 

0 

0 

28 

155 

113 


The number of expected events for each search region in the multilepton channel is given in Table 12 
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Table 12: LHC Multilepton search: Expected events for an integrated luminosity of C = 300 fb -1 . The 
definition of the search channels (ch) is given in the text and Fig. [6] EWK refers here to electroweakly 
produced SUSY particles. 


Ch 

Total SM ± unc. 

STC8 

Zbi 

Xi xt 

X 2 Xm 

Other EWK 

No EWK 

mp+p- < 75 GeV 

1 

10900 ± 3300 

88 

<0.5 

76 

7 

3 

2 

2 

5900 ± 960 

130 

<0.5 

110 

10 

10 

0 

3 

1390 ± 340 

140 

<0.5 

110 

20 

10 

10 

4 

1290 ± 210 

26 

<0.5 

16 

7 

2 

2 

5 

348 ± 121 

19 

<0.5 

10 

6 

2 

1 

6 

45.1 ± 32.4 

8.5 

<0.5 

2.3 

4.5 

0.5 

1.2 

7 

469 ± 125 

26 

<0.5 

9 

12 

4 

2 

8 

29.6 ± 6.8 

9.6 

0.9 

2.5 

5.7 

0.6 

0.8 

9 

1.26 ± 0.41 

1.0 

<0.5 

0.2 

0.6 

0.0 

0.3 

10 

21.4 ± 3.2 

6.6 

1.0 

1.2 

4.6 

0.6 

0.1 

11 

4.48 ± 1.72 

2.7 

0.6 

0.5 

1.9 

0.3 

0.0 

12 

0.0262 ± 0.0095 

0.3 

<0.5 

0.0 

0.15 

0.15 

0.0 

13 

1.06 ± 0.19 

1.1 

0.6 

0.0 

0.4 

0.6 

0.0 

14 

0.89 ± 0.263 

0.3 

<0.5 

0.0 

0.3 

0.0 

0.0 

15 

0.0137 ± 0.0048 

0 

<0.5 

0 

0 

0 

0 


75 GeV < m f+f - < 105 GeV 


1 

111000 ± 16000 

97 

<0.5 

79 

11 

5 

3 

2 

45900 ± 7700 

170 

<0.5 

140 

20 

10 

10 

3 

7490 ± 1390 

210 

<0.5 

140 

50 

10 

10 

4 

4640 ± 490 

26 

<0.5 

12 

10 

2 

1 

5 

994 ± 278 

31 

<0.5 

13 

13 

4 

1 

6 

55.4 ± 40.3 

16 

<0.5 

2 

10 

2 

2 

7 

444 ± 75 

30 

<0.5 

9 

17 

2 

1 

8 

26.2 ± 5.0 

26 

2.9 

6 

16 

4 

2 

9 

1.91 ± 0.47 

2.0 

1.0 

0.2 

1.5 

0.0 

0.4 

10 

16.4 ± 1.4 

6.7 

1.4 

0.9 

4.7 

0.8 

0.4 

11 

5.01 ± 0.93 

3.0 

1.0 

0.6 

1.6 

0.0 

0.8 

12 

0.058 ± 0.0176 

0.12 

<0.5 

0.0 

0.0 

0.0 

0.12 

13 

1.77 ± 0.2 

0.15 

<0.5 

0.0 

0.15 

0.0 

0.0 

14 

2.32 ± 0.35 

0.31 

<0.5 

0.0 

0.31 

0.0 

0.0 

15 

0.113 ± 0.034 

0 

<0.5 

0 

0 

0 

0 

m p+p- > 105GeV 

1 

2380 ± 320 

22 

<0.5 

11 

7 

3 

1 

2 

1720 ± 240 

34 

<0.5 

18 

10 

4 

2 

3 

614 ± 157 

61 

<0.5 

21 

24 

10 

6 

4 

217 ± 47 

10 

<0.5 

3 

6 

1 

1 

5 

57.6 ± 11.5 

11 

0.7 

2 

6 

3 

1 

6 

13.7 ± 3.5 

8.0 

1.2 

1.4 

3.9 

1.8 

0.8 

7 

70.2 ± 6.0 

15 

1.4 

1 

13 

2 

0 

8 

12.6 ± 2.7 

15 

2.6 

2 

10 

2 

1 

9 

0.812 ± 0.174 

1.4 

1.0 

0.2 

0.7 

0.2 

0.4 

10 

11.7 ± 1.4 

3.5 

0.8 

0.2 

2.3 

1.1 

0.0 

11 

1.94 ± 0.41 

2.8 

1.4 

0.3 

1.5 

0.3 

0.7 

12 

0.0391 ± 0.0114 

0 

<0.5 

0 

0 

0 

0 

13 

1.36 ± 0.14 

0 

<0.5 

0 

0 

0 

0 

14 

0.674 ± 0.176 

0.78 

<0.5 

0.0 

0.15 

0.15 

0.47 

15 

0.0137 ± 0.0055 

0.24 

<0.5 

0.0 

0.0 

0.0 

0.24 
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